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-  This  study  is  a  continuation  of  an  engineering  seismology  research  effort 
prompted  by  the  sensitivity  of  guidance  sets  in  Minuteman  Wing  V  to  dis¬ 
tant  earthquakes.  An  earlier  report  considers  the  probable  cause  of 
anomalous  patterns  of  seismic  alarms  triggered  by  two  North  American  earth¬ 
quakes.  This  report  extends  the  previous  study  by  examining  the  propagation 
of  surface  waves  from  the  1975  Pocatello  Valley,  Idaho  earthquake  sequence 
across  Wyoming  to  Wing  V.  In  addition,  the  more  general  question  of 
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surface  wave  phenomena  in  sedimentary  basins  is  addressed,  particularly  the 
effect  of  laterally  inhomogeneous  (dipping)  basin-bedrock  interfaces.  Find¬ 
ings  indicate  that  fundamental  and  first  overtone  surface  waves  are  signifi¬ 
cantly  modified  by  the  travel  path.  In  contrast,  higher  modes  are  relatively 
unchanged  by  the  travel  path,  and  affect  Wing  V  in  much  the  same  way  as  body 
waves  considered  in  the  previous  study. 
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SECTION  1 


INTRODUCTION 

This  work  is  a  continuation  of  an  engineering  seismology  research  effort 
prompted  by  the  observed  sensitivity  of  Minuteman  Wing  V  to  distant  earthquakes. 
An  earlier  report,  "Site  Dependent  Ground  Motions  from  Distant  Earthquakes," 
Wojcik  et  al.  [1980],  considered  the  probable  cause  of  anomalous  patterns  of 
seismic  alarms  triggered  by  two  North  American  earthquakes — 1975  Pocatello 
Valley,  Idaho  (M^  ■  6)  and  1979  St.  Elias,  Alaska  (Mg  *7.1).  It  concluded 
that  body  waves  amplified  by  the  2-3  km  sedimentary  layer  under  Wing  V,  near 
longer  period,  lightly  damped  resonances  of  the  missile  suspension  system  were 
probably  responsible  for  the  patterns.  Layered  half space  calculations  indi¬ 
cated  that  peak  ground  shaking  in  the  period  window  of  Interest  could  vary  by 
a  factor  of  two  at  adjacent  missile  flights. 

The  present  research  report  extends  the  previous  study  by  examining  the 
propagation  of  Love  and  Rayleigh  surface  waves  from  the  Pocatello  Valley, 

Idaho  event  across  Wyoming  to  Wing  V.  In  addition,  the  more  general  question 
of  surface  wave  phenomena  in  sedimentary  basins  is  addressed,  particularly 
the  effect  of  laterally  inhomogeneous  (dipping)  basin-bedrock  interfaces. 

The  report  includes  a  reconnaissance  of  geology,  tectonics  and  velocity  func¬ 
tions  on  the  surface  wave  travel  paths;  an  examination  of  relevant  seismograms 
and  associated  dispersion  calculations  to  pick  the  "best"  crustal  model;  and 
ray  theory  analyses  to  show  the  effects  of  travel  path  inhomogeneities  on 
natural  surface  waves. 

1.1  Travel  Path  Reconnaissance  and  Crustal  Velocity 

The  700  km  travel  path  from  the  Pocatello  Valley  epicenter  to  Wing  V  (70 
km  east  of  Cheyenne,  Wyoming)  includes  180  km  of  Che  Idaho-Wyoming  over thrust 
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belt,  350  km  of  the  Green  River  Basin  of  Wyoming  and  about  170  km  of  the 
Laramie  Range  and  Denver  Basin  (see  Fig.  1-1).  It  is  structurally  complex, 
including  inhomogeneities  in  the  upper  crust  (deep  sediments  over  Irregular 
bedrock,  faults  and  mountainous  uplifts),  as  well  as  variable  depth  to  Moho 
(across  the  Laramie  Range)  and  variable  upper  mantle  velocities  (based  on  Pr 
earthquake  phases) . 

Figure  1-2  shows  travel  paths  from  the  Pocatello  Valley  epicenter  to 
Wing  V,  and  from  the  epicenter  to  the  Golden,  Colorado  seismographic  station. 
Sediment  depth  over  these  two  paths  are  shown  in  Fig.  1-3.  A  number  of  seismic 
surveys  have  been  conducted  near  these  paths.  Most  relevant  are  refraction  sur¬ 
veys  made  by  the  USGS  which  cross  one  or  both  paths  as  shown.  The  western  line 
runs  southeast  from  the  Snake  River  Downwarp  (Fig.  1-1)  across  the  Overthrust 
Belt  to  Flaming  Gorge,  near  the  Utah-Wyoming  border.  The  eastern  survey  runs 
south  from  Sinclair,  Wyoming  through  Wolcott,  Colorado,  paralleling  the  Central 
and  Southern  Rocky  Mountains. 

Velocity  functions  are  shown  in  Fig.  1-4  for  the  Flaming  Gorge  model  in 
western  Wyoming  and  the  Sinclair  and  Wolcott  models  to  the  east.  These  indi¬ 
cate  a  nominal  crustal  thickness  of  40  km  which  thickens  to  about  48  km  across 
a  northeast  trending  shear  zone  indicated  in  Fig.  l-l . 

1.2  Seismograms  and  Theoretical  Interpretation 

Seismograms  of  the  1975  Pocatello  Valley  earthquake  sequence  were  recorded 
on  bedrock  at  Golden,  Colorado  (see  Fig.  1-5).  These  are  the  closest  available 
records  to  Wing  V  and  are  expected  to  suggest  its  ground  motion  history.  The 
main  shock  drove  the  instruments  off  scale  for  over  five  minutes  providing  no 
useful  data,  but  numerous  fore  and  aftershocks  were  well  recorded.  The  short 
period  seismograms  showed  the  usual  sequence  of  arrivals  for  a  local  event, 
e.g.  P  ,  P  ,  S  /L  but  detailed  period  data  could  not  be  inferred.  However, 

n  8  8  8 
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a  long  period  seismogram  of  a  ■  4.7  aftershock,  shown  in  Fig.  1-5,  provided 
useful  dispersion  data.  The  wave  arrival  time  and  particle  motion  polariza¬ 
tion  indicated  chat  these  were  fundamental  mode  Rayleigh  waves. 

From  the  dispersed  train  of  surface  waves  shown  in  Fig.  1-5,  dispersion 
(arrival  time  vs.  period)  and  particle  motion  (ellipticity)  data  were  deduced. 
These  are  shown  in  Fig.  l-6a.  The  retrograde  elliptical  path  identifies  the 
phase  as  a  Rayleigh  wave.  The  measured  group  velocity  data  were  used  to  find 
the  "best”  average  travel  path  model.  The  theoretical  first  and  second  mode 
dispersion  curves  shown  were  for  the  Flaming  Gorge  crust  with  2  km  of  sediment. 
For  comparison,  the  first  mode  for  zero,  two  and  five  kilometers  of  sediment 
is  drawn  in  Fig.  l-6b.  Figure  l-6a  indicates  a  group  velocity  minimum  corres¬ 
ponding  to  Airy  phases  may  be  present  in  the  2  to  4  second  period  range,  to 
which  Ulng  V  has  proven  sensitive. 

1.3  Seismic  Rays  and  Surface  Waves 

Surface  waves  are  usually  analyzed  from  a  modal  viewpoint.  In  laterally 
inhomogeneous  structure  a  ray  theory  perspective  is  more  useful  because  rays 
are  easily  continued  through  Inhomogeneities,  while  modal  solutions  are  not. 

A  suite  of  ray  fans  in  the  Flaming  Gorge  crust.  Fig.  1-4,  are  drawn  in  Figs. 
l-7a,b  for  various  penetration  scales.  Accompanying  these  are  reduced  travel¬ 
time  curves.  Figure  l-7a  shows  mantle  and  crustal  P-rays  while  Fig.  l-7b  shows 
sedimentary  P-rays.  A  correspondence  between  such  rays  (including  S-rays  which 
are  not  shown)  and  surface  wave  phases  exists  and  can  be  used  to  deduce  the 
depth  of  component  ray  penetration  as  a  function  of  surface  wave  period.  The 
finding  is  that  fundamental  mode  Love  and  Rayleigh  waves  and  the  first  Rayleigh 
overtone  in  the  3-4  second  period  range  of  importance  to  Wing  V  are  composed 
of  rays  which  penetrate  to  relatively  shallow  depths,  e.g.  on  the  order  of 
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7  km  or  less.  Therefore,  these  surface  waves  would  be  expected  to  Interact 
with  near  surface  lateral  inhomogeneities  shown  in  Fig.  1-3. 

To  determine  what  effect  near  surface,  lateral  inhomogeneities,  such  as 
the  3“-5#  updip  in  Fig.  l-3b,  have  on  shallow  surface  waves,  two-point  ray 
tracing  calculations  have  been  performed  in  prismatic  structure.  An  example 
is  shown  in  Fig.  l-8a  for  a  layer  over  half space,  modeling  the  updipping  geol¬ 
ogy  in  Fig.  l-3b  as  well  as  a  corresponding  downdip.  The  input  is  a  fundamen¬ 
tal  mode  Love  wave  ilor  a  range  of  periods.  For  updip  the  surface  response 
amplitude  as  a  function  of  distance  over  the  transition  shows  amplification 
followed  by  a  rapid  cutoff.  Amplification  is  due  to  the  converging  geometry 
and  the  cutoff  occurs  because  the  component  rays  leak  into  the  lower  half space, 
to  emerge  elsewhere  depending  on  their  takeoff  angle.  A  similar  suite  of  plots 
is  drawn  in  Fig.  l-8b  for  different  dip  angles  and  waveguide  thicknesses,  for 
comparison.  The  qualitative  features  are  similar. 

The  most  pronounced  effect  of  Love  waves  into  an  updipping  transition  is 
the  distinct  cutoff  at  some  distance,  depending  on  surface  wave  period.  The 
cutoff  is  due  to  energy  leakage  into  the  lower  halfspace  and  the  question  is 
where  does  this  energy  go.  This  is  useful  to  know  because  it  involves  the  idea 
of  detection — particularly  the  question  of  deducing  activity  within  a  basin 
from  seismic  signals  recorded  outside.  Figure  1-9  shows  the  mechanism  of  leak¬ 
age.  For  Love  waves  into  an  updip,  the  energy  is  tramsmltted  as  distinct  ray 
beams  into  the  halfspace,  By  ray  tracing,  it  is  simple  to  determine  where 
these  beams  will  re-emerge  on  the  surface. 

The  results  shown  above  for  Love  waves  propagating  into  prisms  in  Fig. 
l-8a,b  are  an  intermediate  result  of  a  more  general  study  addressing  the  ques¬ 
tion  of  Rayleigh  and  Love  type  ducted  surface  waves  within  a  basin  Interacting 
with  dipping  bedrock  interfaces.  Such  interfaces  are  common  near  basin  edges. 
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as  well  as  under  smooch  depch  transitions  In  the  basin  interior.  Fig.  1-10. 

The  fill  is  idealized  as  a  dissipative  homogeneous,  isotropic  elastic  medium. 
Underlying  bedrock  is  similarly  elastic,  but  with  higher  seismic  wavespeeds. 
Scale  of  the  basins  may  be  on  the  order  of  tens  to  hundreds  of  kilometers. 

This  idealization  provides  a  prototype  model  from  which  fundamental  sur¬ 
face  wave  reflection  and  transmission  phenomena  can  be  recognized  and  quanti¬ 
fied.  Relevant  effects  include:  patterns  of  amplification  and  attenuation 
over  the  free  surface;  energy  leakage  and  its  dependence  on  interface  geometry 
and  surface  wave  type;  influence  of  material  damping,  fill  depth,  impedance 
contrast,  etc.  The  purpose  of  the  analysis  Is  to  provide  both  qualitative  and 
quantitative  data  for  the  analyst  who  must  assess  anomalous  wave  effects  and 
seismic  hazard  due  to  surface  waves  in  real  geologies. 


1 .A  Prototype  Problems  and  the  Seismic  Prism  Analysis  Method 

Three  prototype  problems  are  examined.  They  are  illustrated  in  Fig.  1-11, 
each  consisting  of  an  upstream  waveguide  supporting  a  single  mode  Love  wave, 
connected  to  a  converging  or  diverging  layer.  Figure  1-1 la  shows  an  updlpping 
layer  which  terminates  at  the  vertex;  Fig.  1-llb  shows  a  downdipping  layer 
which  does  not  terminate;  and  Fig.  1-llc  illustrates  the  more  general  case  of 


a  dipping  layer  connecting  waveguides  of  different  depth,  which  includes  the 
ocher  two  as  special  cases  (and  is  not  analyzed) .  These  models  exhibit  most 
relevent  physics  of  the  natural  features  at  some  level  of  idealization. 


In  general,  prototype  models  are  very  useful  in  seismic  analysis  of  complex 
structure.  For  example,  Love's  prototype  waveguide,  so  important  in  early  seis¬ 
mology,  was  the  homogeneous  layer  over  a  half space.  Love  [1911].  The  examples 
examined  here  are  the  simplest  generalization  of  this  basic  model  to  laterally 
inhomogeneous  waveguides. 
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The  mathematical  analysis  of  surface  waves  in  the  prototype  models  is 
based  on  a  ray  decomposition  of  the  incident  surface  wave,  yielding  a  two-point 
ray  tracing  problem  in  a  so-called  seismic  prism.  The  prism  models  the  dipping 
layer  and  is  analyzed  on  the  basis  of  self-similarity.  An  exact  solution  is 
found  for  all  Love-wave-derived  rays  using  a  canonical  solution  for  rays  in 
the  prism.  The  Rayleigh  wave  ray  solution,  which  may  involve  thousands  of  P- 
and  SV-waves,  is  conveniently  displayed  by  a  binary  tree  data  structure.  The 
complete  ray  theory  has  yet  to  be  implemented. 

Secondary  rays,  associated  with  diffraction,  are  free  surface  Rayleigh 
waves,  head  waves  and  cylindrical  waves.  These  are  not  included  in  asymptotic 
ray  theory.  To  include  them  requires  a  canonical  solution  for  edge  diffraction 
in  elastic  media.  The  continuation  of  such  diffracted  solutions  away  from  the 
diffractor  proceeds  by  ray  tracing,  as  in  the  geometrical  theory  of  diffraction, 
Keller  [1958]. 

1.5  Summary 

There  are  three  basic  questions  motivating  the  research: 

1.  What  was  the  causative  wave  type  for  Wing  V  alarms? 

2.  What  are  the  response  anomalies  due  to  surface  waves  through 
inhomogeneous  waveguides? 

3.  How  does  surface  wave  energy  leak  out  of  basins  and  where  does 
it  go? 

Question  1  relates  directly  to  the  Wing  V  alarm  problems  from  the  Poca¬ 
tello  Valley  quake.  Travel  path  from  Pocatello  Valley  to  Wing  V  was  through 
a  number  of  geologies  across  Wyoming,  all  characterized  by  thinning  of  sedimen¬ 
tary  cover  in  the  direction  of  propagation.  Gross  features  described  in  a 
previous  report  indicate  convergence  in  both  azimuth  and  depth  (funneling) . 
Previous  work  does  not  explain  what  effects  this  has  on  surface  waves  on  the 


period  range  of  interest  (1-5  sec).  The  present  work  indicates  that  the  updip 
across  Wyoming  may  filter  out  certain  1-5  second  waves. 

Question  2  addresses  the  problem  of  anomalous  response  in  a  nonuniform 
basin  due  to  surface  waves.  For  example,  the  effects  on  surface  waves  at 
Wing  V  of  the  highland/lowland  depth  change  and  of  interface  steps  or  trans- 
tlons  are  examined  below  as  well  as  the  effects  of  the  Western  Rocky  Mountain 
flank  on  surface  waves.  The  results  bear  on  the  general  question  of  basing 
sensitive  structures  near  basin  edges  or  near  other  waveguide  inhoaogeneities. 

Question  3  concerns  the  seismic  detection  problem.  For  a  source  within 
a  basin  (moving  or  stationary)  most  energy  is  trapped  except  for  direct  body 
waves  which  escape.  We  find  that  the  trapped  surface  wave  energy  is  projected 
out  of  the  basin  as  diffraction  limited  ray  beams  at  the  updipping  edges  or 
transitions.  In  the  sequel  these  effects  are  quantified  on  the  basis  of 
geometry,  wave  type,  impedance  contrasts,  etc.  Typical  take-off  angles  and 
crustal  ray  paths  will  be  described  for  these  leaky  beams. 


Figure  1-1.  Tectonic  nap  of  Wyoming  and  adjacent  states  showing  principal  uplifts  and  basins.  The 
Pocatello  Valley  epicenter  is  marked  by  a  star  on  the  border  above  the  Great  Salt  Lake 
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Figure  l-3a.  Sediment  depth  on  500  km  of  the  travel  path 
from  Pocatello  Valley  to  Wing  V. 
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Figure  l-3b.  Sediment  depth  on  250  km  of  the  travel  path  from 
Pocatello  Valley  to  the  Golden,  Colorado  seismo- 
graphlc  station. 


Vertical  Component 


North- South  Component 

/^\'VWK;W 


East -West  Component 

^'VW\A/wwvww^*vv. 

Pocatello  Valley  aftershock:  origin  time  13h1m20* 


Figure  1-5.  Long  period  WWSSN  seismograms  of  a  Ml  ■  4.7, 
March  29,  1975  Pocatello  Valley  aftershock 
recorded  at  Golden,  Colorado. 
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Figure  l-6a.  Observed  and  theoretical  Rayleigh 
wave  dispersion  data. 
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Figure  l-6b.  A  comparison  of  first  modes  for  0,  2  and  5  km 
of  sediment  over  the  Flaming  Gorge  crustal 
model . . 


Figure  l-7b.  Shallow  penetrating  sediment  ’’-rays  in  the  Flaming  Gorge  crust  with  3-5  km  of  sediments 


Figure  l-8fa.  Fundamental  mode  surface  response  for  the  5  km  layer 
and  various  dip  angles. 


a.  Converging,  terminating  waveguide 


b.  Diverging,  non-terminating  waveguide 


c.  Transitional  waveguide 


Figure  1-11.  Prototype  waveguide  models  to  be  examined 
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SECTION  2 


GEOLOGICAL  AND  SEISMOLOGICAL  RECONNAISSANCE 

This  section  presents  a  reconnaissance  of  geologic,  tectonic  and  seismic 
data  in  the  neighborhood  of  surface  wave  travel  paths  from  the  Pocatello  Valley 
epicenter  to  both  Wing  V  and  a  seismographic  station  near  Golden,  Colorado.  It 
includes  a  description  of  near-surface  structure,  lateral  inhomogeneities  in  the 
deep  sediment-bedrock  interface,  and  velocity  models  for  sedimentary  and  crustal 
sections.  Seismograms  of  the  Pocatello  Valley  event,  recorded  on  Rocky  Mountain 
bedrock  150  km  south  of  Wing  V,  are  described  and  reduced  to  crustal  dispersion 
data  for  fundamental  mode  Rayleigh  waves  with  periods  greater  than  5  seconds.  In 
conjunction  with  theoretical  dispersion  calculations,  the  seismographic  data  are 
used  to  deduce  dispersion  curves  in  the  period  window  of  interest— 1  to  5  seconds. 

2.1  The  Regional  and  Seismic  Setting 

Wing  V  is  distributed  over  the  tri-state  area  of  Wyoming,  Nebraska  and  Colo¬ 
rado  in  the  northern  Denver  Basin  (Julesberg  Basin)  between  the  North  and  South 
Platte  Rivers.  It  is  bounded  to  the  west  by  the  Laramie  and  Front  Ranges  of  the 
Central  and  Southern  Rocky  Mountains,  and  to  the  east  by  the  Great  Plains  of 
Central  North  America.  The  areal  extent  and  neighboring  tectonic  features  are 
illustrated  in  Fig.  1-1,  King  [1969],  showing  principal  uplifts  and  sedimentary 
basins.  Structure  to  the  west  is  dominated  by  the  Cordllleran  system  of  mow- 
tains,  basins  and  plateaus  extending  800-1600  kilometers  inland  from  the  Pacific 
Coast  along  the  length  of  North  America.  Of  principal  interest  here  is  the 
eastern  part  of  the  system  which  was  deformed  most  recently  during  the  Laramide 
orogeny  towards  the  end  of  the  Mesozoic  era.  King  [1977]. 

The  Pocatello  Valley  earthquake  sequence  of  late  March  and  April,  1975  was 
centered  around  the  main  shock  (M^  *6.0)  in  an  area  just  north  of  the  Idaho- 


Utah  border  at  42.1*N,  112. 5°W.  The  locus  of  epicenters,  Arabasz  et  al.  [1981] 
is  indicated  by  a  star  above  the  Great  Salt  Lake  in  Fig.  1*1.  Historically, 
this  is  a  region  of  high  seismicity  within  the  southern  Intermountain  Seismic 
Belt,  Smith  and  Sbar  [1974].  The  most  significant  recent  events  in  the  area  were 
the  1934  Hansel  Valley  (M  -  6.6)  and  the  1962  Cache  Valley  (M  ■  5.7)  earthquakes, 
both  within  75  km  of  the  Pocatello  Valley  epicenter.  Structurally,  the  Poca¬ 
tello  Valley  area  is  a  10  x  15  1cm  northern  trending  sediment  filled  graben  on 
the  eastern  margin  of  the  Basin  and  Range  province.  Arabasz  et  al.  [1981]  con¬ 
clude  that  the  1975  earthquake  sequence  accompanied  irregular  graben  subsidence 
within  the  valley. 

2.2  Travel  Path  Geology  and  Tectonics 

Travel  paths  from  the  Pocatello  Valley  epicenter  (POC)  to  Wing  V  (WG5) , 

and  to  the  Golden  seismographic  station  (GOL)  are  drawn  in  Fig.  1-2.  The  700  km 

/ 

path  from  POC  to  WG5,  across  the  lower  third  of  Wyoming,  includes  about  180  km 
of  the  Idaho- Wyomin g  Over thrust  Belt,  350  km  of  the  Green  River  Basin  (across 
the  Continental  Divide)  and  170  km  of  the  Laramie  Range  and  Denver  Basin  (see 
Fig.  1-1).  The  650  km  path  from  POC  to  GOL  includes  perhaps  150  km  of  the 
southerly  extension  of  the  Over thrust  Belt,  350  km  of  the  Green  River  Basin 
paralleling  the  Uinta  Mountains  and  150  km  across  the  Front  Range. 

The  mountains  and  uplifts  near  the  travel  paths  in  Wyoming  and  Colorado 
were  affected  most  recently  by  the  Laramide  orogeny,  beginning  in  late  Creta¬ 
ceous  time.  The  orogeny  created  great  vertical  uplifts  and  the  present  struc¬ 
tures  developed  from  rocks  of  the  ancient  continental  platform.  Heavily  meta¬ 
morphosed,  steeply  tilted  precambrian  rocks  form  most  of  the  mountain  cores.  A 
discontinuity  in  rock  type  occurs  across  the  northeast  trending  shear  zone 
(Jtollen  Creek-Nash  Fork  Shear  Zone)  shown  in  Fig.  1-2.  The  overall  surface 
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structure  has  since  been  modified  by  processes  during  Tertiary  time.  These 
Include  basin  formation  and  regional  uplift  responsible  for  the  present  height 
of  the  Rockies.  Streams  invigorated  by  the  uplift  have  etched  out  the  ranges 
and  filled  the  basins. 

The  depth  of  basin  sediments  (late  Paleozoic  to  early  Tertiary)  covering 
Precambrian  basement  rock  is  shown  in  Pig.  l-3a,b  for  segments  on  both  travel 
paths  Indicated  in  Fig.  1-2.  These  are  drawn  from  state  topographic  maps  and 
the  tectonic  map  of  North  America.  The  cross  sections  show  non-uniform  cover, 
generally  chinning  towards  the  east  with  the  interface  dipping  intermittently 
on  the  order  of  2°  to  10s.  The  P0C-WG5  path  is  the  more  variable  of  the  two, 
exhibiting  both  updips  and  downdlps  as  well  as  significant  azimuthal  variations 
The  variability  with  azimuth  forms  a  crude  three-dimensional  channel  bounded 
abruptly  by  the  Wind  River  Range  and  associated  faults  and  uplifts  to  the  north 
converging  into  the  Hanna  Basin.  The  POC-GOL  path  is  more  uniform  but  with  an 
80  km  5°  updip  terminating  on  Che  western  flank  of  the  Rocky  Mountains.  The 
most  significant  azimuthal  variations  occur  on  the  flank  of  the  Uinta  uplift  to 
the  south. 

2.3  Seismic  Surveys  and  Velocity  Models 

Seismic  velocity  data  near  the  P0C-WG5  and  POC-GOL  travel  paths  are  avail¬ 
able  from  a  variety  of  sources  including  oil  well  sonic  logs  (near  surface)  and 
seismic  reflection  and  refraction  surveys  (deeper  structure).  Refraction  and 
reflection  surveys  conducted  by  USGS  and  C0C0RP,  respectively,  are  indicated  in 
Fig.  1-2.  Of  principal  use  here  are  reinterpretations  of  the  USGS  data  by 
Prodehl  [1979]  and  Prodahl  and  Pakiser  [1980]  using  smooth  velocity  functions 
rather  than  the  conventional  but  approximate  piecewise  constant  functions.  The 
western  refraction  line  in  Fig.  1-2  runs  southeast  from  the  Snake  River  Down- 
warp,  across  the  Overthrust  Belt  to  Flaming  Gorge  near  the  Utah-Wyoming  border. 
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The  eastern  refraction  line  runs  south  from  Sinclair,  Wyoming  through  Wolcott, 
Colorado,  paralleling  the  Central  and  Southern  Rocky  Mountains.  Of  more  quali¬ 
tative  use  are  the  COCORP  deep  reflection  profiles  acros.  the  Wind  River  and 
Laramie  Ranges  in  Wyoming.  These  provide  a  striking  picture  of  tilted  pris¬ 
matic  structure  in  the  Ranges  and  surrounding  deep  sediment,  e.g.  Smithson  et 
al.  [1979] ,  as  well  as  mapping  a  rapid  increase  in  Moho  depth  across  the  Hash 
Fork-Mullen  Creek  Shear  Zone,  Allmendinger  [1980] . 

Velocity  functions  are  shown  in  Fig.  1-4  near  Flaming  Gorge  on  the  western 
line,  and  Sinclair  and  Wolcott  on  the  eastern.  The  cross  hatching  indicates 
questionable  values  associated  with  possible  low  velocity  zones  in  the  crust, 
and  the  mantle  velocity  gradient  below  the  Moho.  The  Flaming  Gorge  and  Sinclair 
models  show  a  depth  to  Moho  of  about  40  km  while  the  Wolcott  Moho  is  near  48  km. 

The  eastern  models  indicate  a  low  velocity  zone  in  the  upper  crust  common  to  the 
region  Immediately  west  of  the  Front  Range  in  Colorado,  Prodehl  and  Pakiser  [1980], 
but  absent  in  the  more  westerly  Flaming  Gorge  model,  Prodehl  [1979].  Upper  crust 
velocities  are  significantly  higher  on  the  eastern  line  in  contrast  to  the  wes¬ 
tern  because  the  latter  survey  sampled  deep  sediments  between  Bear  Lake  and 
Flaming  Gorge,  while  the  former  sampled  Rocky  Mountain  bedrock  between  Sinclair 
and  Wolcott.  Note  in  Fig.  1-2  that  the  eastern  line  crosses  the  northeast  trend¬ 
ing  shear  zone  which  marks  a  discontinuity  in  surface  rock  type. 

Although  the  Wolcott  velocity  function  south  of  the  shear  zone  exhibits  a 
more  complicated  low  velocity  zone  compared  to  Sinclair,  as  well  as  higher 
near-surface  velocity,  the  principal  difference  is  depth  to  Moho.  Based  on  the 
Laramie  Range  COCORP  reflection  line  data,  Allmendinger  [1980]  found  that  the 
depth  transition  correlates  with  the  Mullen  Creek-Hash  Fork  Shear  Zone  indicated. 
This  implies  that  the  Shear  Zone  is  a  crustal  rather  than  surficial  feature, 
possibly  extending  into  the  mantle.  Therefore,  the  Wolcott  model  is  presumed 
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Co  be  representative  of  crustal  structure  east  of  the  Shear  Zone,  particularly 
under  the  Denver  Basin.  It  is  consistent  with  a  layered  interpretation  of  CJSGS 
refraction  data  in  eastern  Colorado  by  Jackson  et  al.  [1963]. 

The  above  crustal  models  do  not  Include  velocities  in  near-surface  sedi¬ 
mentary  rock.  They  can  be  augmented,  however,  employing  sonic  oil  well  logs  in 
the  shales  and  sandstones  ubiquitous  to  the  basins  on  and  around  the  travel 
paths  in  Fig.  1-2.  The  most  significant  recent  basin  fills  are  of  Eocene  age, 
deposited  after  the  climax  of  the  Laramide  orogeny.  For  example,  the  Green 
River  Basin  contains  thick  Eocene  beds  of  oil  shale.  These  lay  unconformably 
on  older  sedimentary  rocks,  chiefly  Mesozoic  and  late  Paleozoic,  King  [1977]. 
Examples  of  sonic  log  data  reduced  by  Thompson  et  al.  [1981]  from  the  Wind  River 
and  Powder  River  Basins  are  plotted  in  Fig.  2-1.  Velocity  depends  on  depth  and 
rock  type;  for  shales  and  sandstones  the  dependence  has  been  extensively  studied, 
e.g.  Haskell  [1941],  Faust  [1951].  Faust's  velocity  fits  of  well-controlled  shale 
and  sandstone  sections  at  various  depths  are  plotted  in  the  figure  for  compari¬ 
son.  The  velocity  function  is  represented  by  Az1^  ,  where  z  is  depth  and  A 
depends  on  geologic  age.  Faust's  constants,  modified  to  yield  velocities  in  km/ 
sec  for  depth  in  meters,  are  A  (Pennsylvanian)  -  1.1  ,  A  (Mesozoic)  -  1.0  , 

A  (Eocene)  -  .9  and  A  (Post-Eocene)  -  .8  . 


2.4  Golden  Seismograms 


The  Pocatello  Valley  earthquake  sequence  included  a  ■  4.2  foreshock, 
a  6.0  mainshock,  one  4.7  and  two  3.8  aftershocks  and  over  50  lesser  events 


greater  than  3.0,  Arabasz  et  al.  [1981].  Long  and  short  period  WWSSN  seismograms 
of  the  events  were  obtained  from  the  Golden  seismographic  station,  operated  by 
the  Colorado  School  of  Mines,  Majors  [1980],  Glover  [1980].  These  are  the  closest 


available  records  within  the  Denver  Basin,  and  were  recorded  on  bedrock  at  the 


33 


base  of  the  Rocky  Mountains.  Because  WG5  and  GQL  are  only  150  km  apart,  with  no 
significant  difference  in  tectonic  province,  the  GOL  seismograms  were  judged 
adequate  to  infer  bedrock  motion  near  WG5. 

Epicentral  distance  to  GOL  is  about  650  km  and  the  azimuth  is  67°  east  of 
south.  The  main  shock  at  2^3lm6S  UTC,  March  28  drove  the  long  and  short  period 


instruments  off-scale  for  4  to  5  minutes;  however,  the  lesser  events  provided 
good  records.  The  short  period  east-west  seismogram  of  the  first  3.8  magnitude 
aftershock  is  shown  in  Fig.  2-2.  The  sequence,  relative  amplitude  and  timing 
of  arrivals  is  common  to  the  4.2  foreshock  and  the  4.7  aftershock  as  well. 

The  first  arrival,  at  la30S  after  the  origin  time  of  16^15m6s  March  28,  is 
the  Pq  phase.  This  is  a  low  amplitude  refracted  wave  associated  with  the 
crust-mantle  transition  or  Moho.  Following  Pq  is  a  higher  amplitude  Pg  phase 
a  direct  wave  propagated  through  the  granite  crust.  Another  arrival,  P*  ,  assoc i 


ated  with  deeper  layers  in  the  crust  is  present  between  P  and  P  .  The  next 

a*  5 

obvious  arrival  is  Sg  ,  the  direct  shear  wave.  Faster  refracted  phases  (e.g. 

S  )  preceed  this  arrival  but  are  not  obvious  on  the  record.  Love  waves  including 


L  follow  S  but  no  distinctions  can  be  made  and  similarly  for  Rayleigh  waves. 
8  8 

The  corresponding  long  period  seismogram  shows  no  evidence  of  this  event. 


The  only  significant  long  period  record  at  GOL  besides  the  mainshock  is  the 
magnitude  4.7  aftershock  shown  in  Fig.  1-5.  Origin  time  is  13hlm209,  March  29. 
The  long  period  arrival  at  13h5m  begins  with  a  period  near  12  seconds  and  is 
down  to  5  seconds  after  lm30S.  The  vertical  and  east-west  components  have  the 
higher  amplitude.  This  behavior,  the  normal  dispersion  and  particle  motion, 
identifies  the  phase  as  a  Rayleigh  wave.  Group  velocity  of  the  12  second  period 
is  about  2.85  km/sec,  which  is  within  the  expected  scatter  for  fundamental  mode 
continental  Rayleigh  waves,  Oliver  [1962].  As  expected,  the  smaller  shocks  do 
not  excite  surface  waves  very  effectively;  however,  the  4.7  aftershock  is  seen 
to  be  a  marginally  effective  source  and  the  mainshock  is  a  very  rich  source. 


2.5  Observed  and  Theoretical  Dispersion  Data 

The  dispersed  train  of  Rayleigh  waves  observed  from  the  magnitude  4.7  after¬ 
shock  is  used  here  to  deduce  dispersion  (group  velocity  vs.  period)  and  particle 
motion  (ellipticity)  data.  These  data  are  compared  to  theoretical  calculations 
for  each  of  the  crustal  models  described  above,  assuming  a  range  of  sediment  thick¬ 
nesses. 

Direct  measurements  of  predominant  period  and  time  of  arrival  off  Fig.  1-5 
yield  group  velocity  versus  period,  plotted  in  Fig.  2-3.  To  compare  this  travel 
path  dispersion  with  continental  and  sedimentary  dispersion  curves  observed  else¬ 
where  in  the  world,  the  scatter  band  given  by  Oliver  [1962]  is  indicated. 

✓ 

Clearly,  despite  structural  complexity  on  the  travel  path,  no  gross  aberrations 
are  observed  for  dispersed  periods  greater  than  5  seconds. 

Particle  motion  orbits,  i.e.  east-west  versus  vertical  motion,  were  also 
measured.  These  are  shown  in  Fig.  2-3  after  correcting  for  azimuth  to  give  radial 
versus  vertical  motion.  Approximate  ellipticities  are  .41  near  6  seconds  and  .69 
near  9  seconds.  Note  that  in  comparison  to  the  vertical  seismogram  in  Fig.  1-5  , 
the  east-west  amplitude  is  seen  to  smoothly  decrease  with  time  and  predominant 
period.  Actually,  the  smooth  decrease  is  due  to  long  period  WWSSN  instrument 
response,  for  which  gain  drops  linearly  from  3000  to  2000  as  period  decreases 
from  9  to  5  seconds.  Therefore,  observed  east-west  (or  radial)  motion  is  nearly 
uniform  over  the  period  range  while  vertical  motion  fluctuates  principally  as  a 
function  of  ellipticity. 

The  period  window  of  Interest  here  is  1  to  5  seconds.  Using  the  observed 
dispersion  data  available  for  periods  greater  than  5  seconds  as  a  control,  dis¬ 
persion  calculations  in  the  crustal  models  topped  by  sedimentary  sections  were 
extended  into  the  shorter  period  range.  For  example,  assuming  sedimentary  layers 
of  Paleozoic  shales  and  sandstones.  Fig.  2-1,  over  the  Flaming  Gorge  crustal 
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model  yields  the  fundamental  mode  Rayleigh  wave  dispersion  curves  in  Fig.  l-6b. 

A  2  km  layer  gives  the  best  fit  to  data,  but  of  course  optimum  depth  depends  on 
sediment  velocity  which  was  chosen  here  only  to  provide  an  upper  bound.  Velo¬ 
cities  could  be  as  much  as  25%  lower  in  post-Eocene  sediments,  as  indicated  in 
Fig.  2-1.  Nonetheless,  these  calculated  dispersion  curves  show  that  sedimentary 
cover  shifts  the  dispersion  curve  to  the  right  and  introduces  a  pronounced  group 
velocity  minimum  in  the  period  range  of  interest. 

i 

The  dispersion  calculations  were  made  assuming  a  ratio  of  P  to  S  seismic 
velocities  of  2  in  the  sediment  and  1.73  in  the  crust,  corresponding  to  commonly 
assumed  Poisson  ratios  of  1/3  and  1/4,  respectively.  The  sedimentary  ratio  is 
suggested  by  laboratory  measurements  in  shales  from  the  Williston  Basin,  Tosaya 
[1981]  and  the  crustal  ratio  is  consistent  with  rock  data,  Press  [1966]  and  sur¬ 
face  wave  velocity  inversion  results  of  Keller,  et  al.  [1976].  Densities  are 
found  from  data  in  Gardner,  et  al.  t 1974 )  for  sedimentary  rock  and  the  Nafe  and 
Drake  empirical  curves  in  Press  [1966];  although  gravity  modeling  data  in  the 
Green  River  Basin,  Smithson,  et  al.  [1979]  indicates  5-10%  lower  densities  in 
basement  rock.  The  computer  code  used  for  the  dispersion  calculations  assumed 
a  piecewise  constant  layered  half space,  Herrmann  [1978]  with  20-30  layers  to 
adequately  resolve  the  velocity  models. 

A  comparison  of  fundamental  and  first  higher  mode  Rayleigh  wave  dispersion 
is  shown  in  Fig.  2-4  for  the  Flaming  Gorge  and  Sinclair  models  covered  by  a 
3  km  Paleozoic  layer.  Group  velocity  minima  for  the  fundamental  are  near  3.5 
seconds  in  both  models,  and  near  4.3  seconds  for  the  first  overtone.  The  pres¬ 
ence  of  these  minima  in  the  period  range  corresponding  to  Wing  V  suspension  sys¬ 
tem  resonances  is  very  suggestive.  Group  velocity  minima  yield  the  Airy  phase — 
associated  with  long  duration  ringing  due  to  late  arrivals  propagating  near  the 
minimum  group  velocity.  Such  long  duration  ringing  near  critical  periods  of  a 
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surface  structure  nay  excite  significant  response  even  if  the  amplitude  of 
ground  shaking  is  relatively  low. 

The  above  results  are  for  an  average  Moho  depth  of  40  km.  Because  depth 
increases  to  48  km  across  the  Shear  Zone,  about  25-30%  of  the  travel  path  is 
through  this  thicker  crust.  Dispersion  curves  for  the  Wolcott  model  are  there¬ 
fore  shown  in  Fig.  2-5  for  no  sediment,  corresponding  to  the  Golden  seismograph 
location,  and  about  2.5  km  of  sediment  suggestive  of  Denver  Basin  fill  west  of 

Golden  and  similar  to  that  under  Wing  V.  These  show  similar  group  velocity 
minima,  implying  that  the  change  in  crustal  thickness  does  not  radically  alter 
the  dispersion.  A  favorable  comparison  of  observed  and  theoretical  ellipticity 
is  also  shown  in  the  figure  for  9  second  waves.  The  calculated  ellipticity  is 
nearly  constant  (-. 73)  for  periods  from  2  to  12  seconds  so  the  observed  ellip¬ 
ticity  at  6  seconds  shown  in  Fig.  2*3  does  not  compare  as  well. 

2.6  Response  of  Denver  Basin  Sediments  to  Love  and  Rayleigh  Waves 

In  order  to  convert  estimates  of  bedrock  motion  to  response  at  the  surface 
of  deeper  sediments,  surface  wave  response  ratios  were  calculated  for  two  sites 
in  the  Denver  Basin.  The  responses  at  Flights  H  and  F  (1800  m,  3400  m  depth  of 
sediment,  respectively,  over  the  Wolcott  model)  of  Wing  V  were  normalized  by 
surface  wave  response  in  the  Wolcott  model.  This  allows  us  to  Investigate  the 
periods  which  are  most  strongly  amplified  by  the  sediments  in  the  Denver  Basin. 
Flights  H  and  P  are  chosen  because  they  bound  the  sediment  depths  for  all  of 
Wing  V. 

We  assume  that  the  surface  vave  energy  flux  is  the  same  on  a  path  from  the 
Rocky  Mountains  to  the  Denver  Basin  and  we  scale  the  mode  shapes  obtained  from 
dispersion  calculations  on  the  basis  of  normalized  total  energy.  If  Iq  is  the 
normalized  energy  in  the  surface  wave  mode  and  U  is  its  group  velocity,  then 

o 
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Che  ratio  of  surface  displacements  at  two  different  locations  is  given 

by 


V°i 


The  assumption  of  uniform  energy  flux  neglects  scattering  ioto  body  waves  and 

other  surface  wave  modes,  which  can  be  included  only  by  resorting  to  numerical 

% 

methods.  Nonetheless,  the  present  results  are  qualitatively  correct  for  the 
cases  considered. 

Figures  2-6a  and  b  show  results  for  radial  and  vertical  components  of  Ray¬ 
leigh  waves,  respectively.  Nominal  amplification  is  on  the  order  of  3  relative 
to  bedrock  response.  Peak  amplification  is  about  4  and  it  occurs  at  a  period  of 
about  2.7  sec  for  Flight  H  (shallow  sediment)  and  at  about  4.4  sec  for  Flight  P 
(deeper  sediment).  In  Fig.  2-7,  a  nominal  amplification  of  4  is  shown  for  the 
Love  wave  case.  Peak  amplification  is  about  5.3  and  occurs  at  a  period  of  about 
4.4  sec.  The  significant  findings  are:  For  Rayleigh  waves  the  amplification 
peaks  occur  within  the  period  window  of  Interest  (1-5  sec)  and,  in  fact,  bracket 
the  critical  oeriod  of  3.5  sec,  Wojcik  et  al.  [1980],  with  a  relative  ratio  of 
peak  to  nominal  amplitude  of  about  1.33;  and  for  Love  waves  the  peaks  occur  near 
the  upper  end  of  the  period  window  with  a  peak  to  nominal  ratio  near  1.25.  In 
both  cases,  the  peak  relative  amplification  decreases,  and  its  period  increases 
as  the  sediment  depth  is  increased  from  the  2  km  to  3  km  sites.  The  fact  that 
peak  amplification  in  the  Rayleigh  wave  case  brackets  the  critical  period  (3.5 
sec)  and  in  the  Love  wave  case  occurs  at  longer  period  suggests  that  Rayleigh 
waves  are  more  likely  to  produce  patterns  of  response  anomalies  chan  are  Love 
waves . 
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Figure  2-4.  A  comparison  of  Rayleigh  wave  dispersion  for 
Sinclair  and  Flaming  Gorge  models  with  3  km 
of  sediment. 
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Figure  2-5.  Dispersion  curves  for  Wolcott  model. 
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Figure  2-6a.  Radial  displacement  ratios  for  fundamental  mode  Rayleigh  waves. 


Figure  2-6b.  Vertical  displacement  ratios  for  fundamental  mode  Rayleigh 
waves  (see  text  for  explanation). 
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DieplacemenC  ratios  for  fundamental  mode  Love  waves  (see  text 
for  explanation) . 


SECTION  3 


SURFACE  WAVES;  A  SEISMIC  RAY  THEORY  PERSPECTIVE 

Shorter  period  surface  waves  with  periods  of  1  to  5  seconds  are  expected 
to  interact  significantly  with  the  laterally  inhomogeneous  near-surface  struc¬ 
ture  on  travel  paths  crossing  the  Green  River  Basin  and  the  Central  Rockies 
described  above.  The  question  is:  How,  and  to  what  degree,  are  these  surface 
waves  filtered  or  otherwise  modified  by  the  travel  path  before  reaching  the 
Denver  Basin?  To  investigate  this  question,  it  is  necessary  to  first  recognize 
that  there  are  two  distinct  types  of  surface  waves  to  be  considered.  First  is 
the  fundamental  mode  Rayleigh  wave  which  depends  only  on  the  presence  of  a 
free  surface  for  its  existence,  decays  exponentially  with  depth  and  is  composed 
of  constructively  interfering  evanescent  F-  and  SV-  waves.  Second  are  the 
Rayleigh  overtones  as  well  as  the  fundamental  Love  mode  and  its  overtones, 
which  are  all  channeled  waves  composed  in  part  or  fully  of  trapped  plane  waves, 
typically  in  a  layer  from  the  free  surface  to  some  depth  of  maximum  wave  pene¬ 
tration. 

In  an  effort  to  explore  the  question  of  surface  wave  interaction  with  basin 
and  mountain  structure,  this  section  examines  their  effects  on  channeled  surface 
waves  in  the  context  of  seismic  ray  theory.  This  is  a  useful  perspective  because 
rays  are  easily  continued  through  inhomogeneities,  while  evanescent  modal  solu¬ 
tions  are  not.  It  is  helpful  to  consider  first  the  fundamental  mathematical 
analysis  of  the  governing  wave  equations  in  a  half space  with  seismic  velocities 
varying  smoothly  in  the  vertical  direction  only.  This  analysis  is  based  on  the 
theory  of  partial  differential  equations  of  mixed  type  (hyperbolic-elliptic)  and 
relies  on  the  method  of  characteristics  for  the  hyperbolic  solution. 
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3.1  Mixed  (Hyperbolic -Elliptic)  Wave  Fields  Composing  Surface  Waves 

Consider  time  harmonic  Crains  of  Love  and  Rayleigh  surface  waves  propagating 
over  an  elastic  halfspace  in  the  positive  x  direction  with  phase  velocity,  V  , 
Fig.  3-1.  The  half space  is  a  one-dimensional  medium  in  that  wave  speeds  are  a 
function  of  depth  only,  a(z)  for  P-waves  and  0(z)  for  SV  and  SH  waves. 
Particle  motion  for  Love  waves  is  antiplane — in  the  y  direction,  while  for 
Rayleigh  waves  motion  is  inplane — in  the  x  and  z  directions.  Normalized 
displacements  on  the  free  surface  may  be  written  as 

Love  waves:  v  (x,t)  *  cos  u)(t-x/V)  (3-1) 

o 


Rayleigh  waves: 


l 


uQ(x,t) 


e  cos  m(t-x/V) 
o 


(3-2) 


|wo(x,t)  »  sin  iii(t-x/V) 

where  u  ,  v  and  w  are  x  ,  y  and  z  surface  displacements,  respectively, 
o  o  o 

w  is  circular  frequency  and  eQ  is  surface  ellipticity  (ratio  of  horizontal  to 
vertical  amplitude). 

The  governing  wave  equations  describing  any  two-dimensional  motion  of  the 
half space  are 


,  .  32v  32v  1  32v 

Antiplane: 


(3-3) 


32<d 


1  oc 

Inplane:  ^  — 


.  -1  Hi 

v  v  (RzY  3t2 


(3-4) 


where  $  and  ij/  are  the  Lame  displacement  potentials  for  which  u  * 
and  w  ■  -|i  +  |i  .  In  order  to  solve  for  the  surface  wave  field  in  the  half¬ 
space,  it  is  necessary  to  reduce  the  number  of  independent  variables  in  these 
partial  differential  equations  (PDE's).  The  usual  approach  is  to  express  them 
as  a  number  of  ordinary  differencial  equations,  via  separation  of  variables  or 
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integral  transforms.  A  useful  alternative,  however,  is  to  attack  the  wave 
equations  directly  by  means  of  wavefront  and  ray  analysis,  described  next. 

Surface  waves  in  Fig.  3-1  are  clearly  stationary  in  a  coordinate  system 
convected  with  phase  velocity,  V  ,  along  the  x  axis.  Accordingly,  if  t  and 
x  are  replaced  through  the  transformation,  t  ■  t  -  px  where  p  =  1/V  ,  then 
for  the  Love  wave  case,  the  governing  equation, (3-3)  becomes 

■55T  *  <ig(2)  ■gxTT  ;  qg  ■  i/B2(z)  -  p2  (3-5) 

which  is  a  PDE  in  only  two  independent  variables.  From  the  general  theory  of 

second-order  partial  differential  equations,  Sommerfeld  [1964],  this  may  be 

characterized  by  type  as  hyperbolic,  parabolic  or  elliptic,  depending  on  the 

sign  of  the  coefficient  qg(z)  •  The  coefficient  changes  sign  at  depth,  hence 

the  PDE  is  of  mixed  type— hyperbolic  when  q«  >  0  ,  elliptic  when  q2  <  0  and 

p  p 

parabolically  degenerate  when  q^  *  0  .  Therefore,  for  a  monotonically  increas¬ 
ing  velocity  function  the  displacement  field  in  the  halfspace  is  divided  into 
two  domains:  hyperbolic  in  a  layer  above  the  so-called  elliptic  horizon  at 

depth  z  ,  where  S(z  )  ■  V  and  q?(z  )  ■  0  ;  and  elliptic  below  this  depth, 
in  m  p  m 

When  low  velocity  zones  exist,  there  may  be  more  than  one  elliptic  horizon 
defining  hyperbolic  channels,  for  example. 

In  the  convected  wave  equation  (3-5)  the  constant,  p  »  1/V  is  the  hori¬ 
zontal  slowness  and  is  commonly  referred  to  as  the  ray  parameter.  The  coeffici¬ 
ent,  qg(z)  is  the  vertical  slowness  and  is  imaginary  below  the  elliptic  horizon. 
This  follows  because,  in  the  definition  of  q|  (see  (3-5)),  the  l/g2(z)  term 
represents  the  relative  contribution  of  inertial  forces  (recall  it  was  the  co¬ 
efficient  of  acceleration  in  the  wave  equation  (3-3)),  while  the  p2  -  1/V2 
term  is  the  relative  contribution  of  elastic  restoring  forces.  Inertial  forces 
dominate  elastic  forces  when  q|  >  0  (hyperbolic) ,  and  conversely  for  elastic 


forces  when  q*  <  0  (elliptic).  Therefore,  only  inertial  wave  fields  (homogene- 
p 

ous  waves)  exist  in  the  hyperbolic  domain  but  not  in  the  elliptic,  where  steadily 
moving  elastic  equilibrium  fields  (inhomogeneous  or  evanescent  waves)  result, 
decaying  exponentially  with  depth. 

For  the  case  of  Rayleigh  waves,  the  wave  equations  (3-4)  in  convected  co¬ 
ordinates  are 

0  *  qa(z)  0  5  **£  *  1/0t2<z>  '  P2 

(3-6) 

“  4U)  0  • 

Assuming  that  both  a(z)  and  6(z)  increase  with  depth,  the  wave  fields  can  be 
divided  into  hyperbolic  and  elliptic  domains  as  before.  Clearly,  both  potentials 
have  their  elliptic  horizons  on  the  free  surface  when  V  <  6(0)  ,  corresponding 
to  fundamental  mode  (or  free  surface)  Rayleigh  waves  (coupled  inhomogeneous  or 
evanescent  wr_vef ields) .  When  6(0)  <  V  <  o(0)  the  i|/  potential  has  its  ellip¬ 
tic  horizon  at  depth,  while  the  <j>  horizon  Is  on  the  surface;  and  when  V  >  a(0) 
both  <j>  and  have  their  horizons  at  depth.  These  latter  cases  correspond  to 
higher  mode  (or  channeled)  Rayleigh  waves  because  at  least  one  hyperbolic  layer 
or  channel  exists,  admitting  inertial  or  homogeneous  surface  waves. 

3.2  Hyperbolic  Wave  Fields  and  Rays 

The  hyperbolic  wave  fields  follow  from  the  convected  wave  equations  (3-3,6) 

using  the  method  of  characteristics.  Recognizing  that  in  the  hyperbolic  domains 

(3-5,6)  are  one-dimensional  wave  equations  in  the  z  direction  with  variable 

slowness  or  velocity,  the  characteristic  equation  is  simply  —  ■  t  q(z)  (where 

q  means  q  or  q  ) .  Integrating  from  the  free  surface  gives 
Cl  D 

z 

T(z)  -  1  /  q(OdS  +  Tq±  (3-7) 


50 


where  tq+  is  the  free  surface  intercept  of  the  characteristic.  These  charac¬ 
teristic  curves  yield  the  wavefronts  in  the  hyperbolic  channel  as  shown  in 
Fig.  3-2.  Through  each  point  of  the  channel,  there  are  two  characteristics  cor¬ 
responding  to  the  ±  signs  in  (3-7) .  At  the  point  of  reflection  on  the  elliptic 
horizon  the  wavefronts  are  either  cusped  (zero  included  angle)  in  the  case  of 
a  smooth  wave  speed  function  (Fig.  3-2a),  or  bent  (finite  angle)  if  the  wave 
speed  is  discontinuous  at  the  horizon  (Fig.  3-2b). 

In  applications,  the  wavefront  trajectories  or  rays  are  a  more  convenient 
means  of  describing  hyperbolic  wave  fields  than  are  characteristics  because 
rays  are  invariant  with  time.  They  follow  from  the  fact  that  a  wavefront  in 
an  Isotropic  medium  is  a  surface  of  constant  phase  moving  normal  to  itself  at 
the  local  wave  speed.  Representing  this  surface  by  the  phase  function,  $(xrz, 
t)  •  0  ,  then  from  the  equation  of  characertistlcs. 


z 

<t(x,z,t)  -  t  J  q(OdC  +  To  -  T  «  ^(x.z)  -  t  (3-8) 

o 

where 

z 

Y(x,  z)  -  p x  i  /  q(Od?  +  To  (3-9) 

o 

is  the  wave  function  or  eiconal.  The  normal  to  a  wavefront  at  some  fixed  time 
is  just  'TV  ,  or 


n  =  STV  •  Pi  i  q(z)k 


(3-10) 


3  3 

where  7  «  i  +  k  gradient  operator.  This  double  valued  vector  is 

the  slowness  (or  refractive  index)  and  makes  an  angle  ±0  with  the  positive 
z  axis  given  by  tan9  •  p/q(z)  which  is  Snell's  law  for  a  one-dimensional 
medium.  Substituting  the  definition  of  q(z)  yields  the  more  familiar 
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form  of  Snell's  law 


sin8 

TtT)  m  p 


(3-11) 


where  c  is  equal  to  a  or  8  for  P  and  SV  or  SH  rays,  respectively. 

Snell's  law  provides  the  relation  between  ray  angle  and  depth,  hence  the 
basis  for  ray  tracing  in  a  one-dimensional  medium.  Referring  to  Fig.  3-3,  along 
the  ray  -  tanS  ,  hence  substituting  tan9  -  p/q(z)  and  integrating. 


x 


P  / 


-ra-.  +  x 

q(5)  o 


(3-12) 


where  xq  Is  the  starting  point  on  the  free  surface.  Travel  time  is  found  from 
-  l/c(z)  with  ds  ■  dz/cos8  ,  whence 


t 


I 


c2(s)q(c>  +  Co 


(3-13) 


and  tQ  is  starting  time.  In  a  hyperbolic  channel,  the  rays  corresponding  to 
the  wavefronts  in  Fig.  3-2  are  trapped  between  the  free  surface  and  elliptic 
horizon  as  shown  in  Fig.  3-4.  The  skip  distance,  X  and  travel  time,  T  ,  on  a 
ray  segment  from  the  free  surface  to  the  elliptic  horizon  and  back  to  the  free 
surface  are 

zm  zm 

*  *  2»  l  •  1 ' 2  l  ■  (3-l4) 

For  a  smooth  velocity  function  (Fig.  3-5a) ,  the  elliptic  horizon  is  a  caustic, 

which  is  tantamount  to  saying  it  is  a  boundary  between  hyperbolic  and  elliptic 
wave  fields.  When  a  ray  is  reflected  at  the  horizon  (touches  the  caustic), 
the  amplitude  reflection  coefficient  is  unity  but  there  is  a  phase  change  of 
v/2  .  For  a  discontinuous  velocity  function  the  reflection  coefficient  and 
phase  change  are  found  from  the  usual  formulas  for  reflection  at  interfaces. 
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3.3  Ray  Tracing  in  Smooth  Crustal  Models 


To  examine  channeled  surface  waves  on  travel  paths  between  Pocatello 
Valley  and  the  Denver  Basin,  it  is  necessary  to  first  trace  ray  paths  In  the 
crustal  models.  This  amounts  to  evaluating  the  Integral  In  (3-12).  The  easi¬ 
est  approach  is  to  approximate  the  crustal  models.  Fig.  1-4,  by  piecewise 
linear  velocity  segments  and  use  the  fact  that  rays  In  a  linear  gradient  are 
simply  arcs  of  circles  (Officer  [1974],  Exercise  6.5) .  More  powerful  algorithms 
for  smoother  velocity  functions  are  available  but  not  warranted  for  the  purposes 
of  this  section.  Mantle  and  crustal  P-rays  in  the  Flaming  Gorge  model  with 
5  km  of  Paleozoic  sediment  are  shown  in  Fig.  l-7a,  including  reduced  travel  time 
curves.  For  reference,  shallow  penetrating  sedimentary  P-rays  in  the  upper 
crust  are  shown  in  Fig.  l-7b  for  smoothly  increasing  sediment  velocity  and 
for  a  more  realistic  sedimentary  column  corresponding  to  Wing  V  in  the  Denver 
Basin. 

The  mantle  ray  fan  in  Fig.  l-7a  starts  at  19.4°  from  the  vertical,  with  a 
.05°  ray  angle  increment,  to  the  final  ray  at  20.4°.  A  velocity  gradient  of 
.005  km/sec/km  is  assumed  in  the  upper  mantle  below  the  Moho  (approximated 
from  Jeffreys  [1929])  as  indicated  in  Fig.  1-4.  This  value  is  supported  approxi¬ 
mately  by  deeper  refraction  studies  near  the  Rocky  Mountains,  Julian  [1970]. 

The  crustal  ray  fan  in  Fig.  l-7a  is  a  continuation  of  the  mantle  fan  from  about 
20*.  at  increments  of  .25*  to  24s.  Whereas  the  mantle  rays  yield  refraction 
like  arrivals  (pseudo  head  waves)  and  a  nearly  linear  reduced  travel  time  curve 
(by  virtue  of  the  assumed  linear  mantle  gradient),  the  crustal  rays  show  much 
more  variety  in  terms  of  penetration  depth  and  travel  times.  The  crustal  re¬ 
duced  travel  time  curve  shows  a  triplication  from  about  150  km  to  300  km,  con¬ 
sistent  with  the  reduced  USGS  refraction  data,  Prodehl  [l979]»  from  which  the 
crustal  velocity  model  was  originally  derived.  Minor  triplication,  for  example 
at  75  km,  is  due  to  the  rapid  velocity  increase  in  the  upper  crust. 
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With  travel  time  curves  available,  two-point  ray  tracing  allows  an  identi¬ 
fication  of  arriving  phases  in  the  Denver  Basin  and  a  correlation  with  depth  of 
ray  penetration.  Consider  rays  connecting  the  Pocatello  Valley  source  and 
Golden,  Colorado,  an  epicentral  distance  of  A  •  650  km  .  Referring  to  travel 
time  curve  edc  in  Fig.  l-7a,  there  is  one  direct  or  single  skip  mantle  ray 
and  a  series  of  multiple  skip  mantle  rays.  Skip  distances  are  Xffi  *  A/m  where 
m  is  the  integral  number  of  skips,  m  -  1,2,  ...  ,  M  with  an  upper  bound  of 

i 

M  ■  4  corresponding  to  X^  *  162.5  km  .  From  the  mantle  travel  time  curve, 

the  direct  ray  arrives  at  90  sec  and  the  last  multiple  at  116  sec.  Respective 

depths  of  penetration  are  60  to  40+  km.  Next,  from  the  crustal  travel  time 

curve  cb  ,  there  is  no  direct  crustal  arrival  but  three  multiple-skip  crustal 

rays  are  present  with  skip  distance,  X  ■  A/m  ,  m  -  2,3,4  .  Note  that  the 

m 

m  »  2  case  corresponding  to  X^  “  325  km  does  exist  although  not  indicated 
in  the  fan  of  Fig.  l-7a.  The  range  of  calculated  arrival  times  are  105  to  116 
seconds  and  the  depths  are  40  to  32  km.  Finally,  from  the  travel  time  curve 
ba  ,  there  are  an  unlimited  number  of  crustal  skips  for  m  -  2,3,... 
starting  from  116  seconds  and  penetrating  to  depths  from  32  km  up  to  the 
free  surface.  For  example,  the  13th  skip  (X^  "  50  km)  arrives  at  145  seconds 
and  penetrates  to  10  km. 

Referring  to  the  Pocatello  Valley  aftershock  recorded  at  Golden,  Fig.  2-2 

shows  the  correlation  of  early  arriving  body  wave  phases  and  segments  of  the 

travel  time  curves  in  Fig.  1-7.  In  particular,  Pn  corresponds  to  direct  and 

multiple  mantle  P-rays  on  curve  edc  ;  P*  to  multiple  deeper  crust  P-rays  on 

curve  be  ;  and  P^  to  multiple  shallower  crust  P-rays  on  curve  ba  .  In  fact, 

the  phase  Identifiers  on  the  seismogram  were  placed  according  to  the  above 

calculated  arrival  time  (P  :  90s,  P*;  105s,  P  116s).  In  contrast  to  the  usual 

n  g 

* 

description  of  P  and  P^  as  direct  waves,  according  to  this  Interpretation 
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they  are  acf-<.Jly  skipping  rays — with  at  least  one  intermediate  reflection  off 
the  free  surface  for  an  epicentral  distance  of  650  km.  Clearly,  the  same 
arguments  can  be  applied  to  travel  time  curves  for  S-rays.  This  is  somewhat 
more  to  the  point  because  S-rays  are  the  principal  constituent  of  higher¬ 
mode  surface  waves.  Unfortunately,  this  is  beyond  the  scope  of  the  present 
investigation. 

Instead  of  examining  S-rays,  attention  is  directed  to  dispersion  curves 
in  Fig.  3-5  showing  phase  velocity  versus  period  for  higher-mode  Rayleigh  waves 
in  the  Flaming  Gorge  crust  with  0,  3  and  5  km  of  overlying  sediment.  For  5  km 
of  sediment  in  the  period  range  1  to  5  seconds,  1.63  km/sec  -  V  -  2.71  km/sec. 
Therefore,  the  elliptic  horizon  for  S-waves  is3km*3m^6  km.  Therefore, 
SV-rays  composing  this  surface  wave  do  not  penetrate  deeper  than  6  km.  If, 
instead,  we  assume  3  km  of  overlying  sediments  we  find  the  elliptic  horizon  for 
S-waves  is2.5km^3m-9  km.  If  we  assume  no  sedimentary  cover,  the  ellip¬ 
tic  horizon  for  S-waves  is  3.5  km  *  3  m  $  21.5  km.  For  3  km  of  sediment,  the 
elliptic  horizon  for  the  second  Rayleigh  wave  overtone  in  the  period  range  1 
to  5  seconds  is  3  km  ^  3  m  -  27  km. 

Specializing  these  results  to  periods  of  3.5  seconds,  which  are  of  parti¬ 
cular  interest  to  Wing  V,  we  find  that  the  elliptic  horizon  is  within  the  sedi¬ 
mentary  layer  for  both  3  km  and  5  km  sedimentary  depths.  However,  with  no 
sedimentary  cover,  rays  penetrate  to  13  km. 

From  this  we  conclude  that  higher  model  Rayleigh  waves  at  periods  of  inter¬ 
est  to  Wing  V  will  interact  with  the  dipping  bedrock  interface  along  the  travel 
path  from  Pocatello  Valley.  In  the  next  section,  we  will  apply  ray  theory  to 
study  this  interaction  in  order  to  identify  the  modes  and  periods  which  can  be 
propagated  along  the  travel  path  to  Wing  V. 
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3.4  Equivalent  Ray  Beams  and  the  Hyperbolic  Solutions 

The  above  results  show  that  surface  waves  are  composed  of  hyperbolic  and 
elliptic  wave  fields.  The  former  can  be  described  by  simple  rays  while  the 
latter  require  diffracted  rays,  Keller  and  Karal  [1964],  and  are  beyond  the 
scope  of  the  present  work.  Hyperbolic  wave  fields  include  SH-rays  for  all 
Love  waves;  SV-rays  for  Rayleigh  waves  when  V  <  ct(0)  ;  and  both  P-  and  SV- 
rays  for  Rayleigh  waves  when  V  >  ct(0)  .  Elliptic  wave  fields,  which  include 
the  free  surface  Rayleigh  wave  case  and  evanescent  P-wave  case  when  V  <  a(0)  , 
cannot  be  described  by  conventional  rays.  They  can  be  treated  by  numerical 
methods,  which  is  the  subject  of  a  current  study  by  the  present  authors.  It 
follows  that  the  present  study  by  conventional  ray  methods  consider  only  a  por¬ 
tion  of  the  surface  wave  field  which  can  occur,  namely,  that  portion  governed 
by  hyperbolic  equations. 

A  convenient  view  of  the  surface  wave  ray  system  in  a  waveguide  is  illus¬ 
trated  in  Pig.  3-6.  Here  a  source  of  monomode  surface  waves  approaches  from 
the  left  and  is  refracted  from  depths  appropriate  to  the  variation  of  wave 
speed  with  depth.  At  intervals,  X  ,  called  the  skip  distance,  the  rays  encoun¬ 
ter  the  free  surface  and  are  reflected  downward  where  the  process  is  repeated. 
However,  when  the  rays  encounter,  instead  of  the  free  surface,  a  block  of 
material  with  the  same  wave  speeds  as  the  surface  layer,  they  are  refracted 
along  straight  rays  as  Fig.  3-6  shows.  By  reciprocity,  we  can  inject  ray 
beams  from  the  block  into  the  half space  where  they  will  propagate  as  surface 
waves,  provided  the  frequency  and  ray  Inclination  satisfy  the  dispersion 
relation.  Figure  3-7  shows  an  example  of  how  this  view  of  surface  waves  can 
be  used  to  determine  the  effect  of  a  lateral  inhomogeneity — in  this  case  a 
dipping  layer  corresponding  to  an  Idealization  of  the  western  flank  of  the 
Rocky  Mountains  on  the  path  between  Pocatello  Valley  to  Golden,  Colorado, 
shown  in  Fig.  l-3b.  This  case  will  be  treated  quantitatively  in  Section  4. 


This  view  of  surface  waves  can  be  used  to  understand  qualitatively  their 
interaction  with  typical  graben  structures  of  Basin  and  Range  valleys  as  well. 
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Figure  3-5.  Phase  velocities  for  first  Rayleigh  overtone  (second 
mode)  in  the  Flaming  Gorge  crust  with  0,  3  and  5  km 
of  sediment. 
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SECTION  4 


SURFACE  WAVES  IN  SEDIMENTARY  PRISMS 

In  this  section  we  address  questions  2  and  3  raised  in  Section  1.5,  as 
follows : 

•  How  does  the  non-uniform  travel  path  between  Pocatello  Valley  and 
the  Denver  Basin  affect  short  period  surface  waves  reaching  Wing  V? 

•  How  does  surface  wave  energy  leak  out  of  basins  and  where  does  it  go? 

The  sedimentary  layer  along  the  travel  path  may  be  idealized  as  a  sequence  of 
sedimentary  prisms  (dipping  layers  over  a  half  space) .  Here  we  analyze  the 
basic  geologic  units  or  prototype  waveguides  as  illustrated  in  Fig.  4-1.  A 
further  idealization  is  that  the  velocity  in  the  sedimentary  prisms  is  constant 
and  that  it  overlies  a  homogeneous  half space  representing  a  granite  basement 
rock. 

Although  we  have  emphasized  Rayleigh  waves  in  previous  sections,  here  we 
concentrate  on  the  Love  wave  case.  This  is  because  both  types  must  eventually 
be  considered  and  the  Love  wave  case  is  simpler  and  embodies  most  of  the  ideas 
which  will  eventually  be  required  to  treat  the  Rayleigh  wave  case.  This  sec¬ 
tion  presents  the  theory  of  ray  tracing  in  a  prism  and  the  results  of  Love  wave 
propagation  through  prototype  models  of  sections  of  the  travel  path.  These 
include  a  range  of  sediment  thicknesses  and  updip  and  downdip  angles.  The  two 
sediment  thicknesses  considered,  5  km  and  2  km,  correspond  to  the  western  and 
eastern  flanks  of  the  Rocky  Mountains,  respectively.  The  Interface  updlps  and 
downdips  of  3*,  5°  and  10°  are  observed  along  the  travel  path. 

We  find  that  updips  tend  to  promote  energy  leakage  into  the  underlying 
halfspace  while  downdips  tend  to  trap  energy.  We  are  able  to  trace  the  leaked 
energy  into  the  crust  and  to  determine  where  it  re-emerges  at  the  surface. 

This  information  is  relevant  to  the  hypothesis  that  the  patterns  of  seismic 
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alarms  at  Wing  V  are  due  Co  leaked  surface  wave  energy  which  is  refracted  under 
the  Rocky  Mountains.  It  is  also  relevant  to  the  general  problem  of  reflection 
and  refraction  of  surface  wave  energy  at  basin  edges.  These  results  help  to 
screen  candidate  sites  for  facilities  which  may  be  sensitive  to  periods  of 
ground  shaking  that  are  conditioned  by  geology  peculiar  to  basin  edges. 

4.1  Two-Point  Ray  Tracing  in  Sedimentary  Prisms 

With  the  surface  wave  decomposed  into  rays,  the  problem  of  calculating 
response  in  a  dipping  region  of  the  waveguide  due  to  first-order  rays  can  be 
reduced  to  ray  tracing.  This  also  involves  phase  and  amplitude  determination, 
followed  by  a  ray  sum  at  the  output  point.  Ray  tracing  is  simplified  consid¬ 
erably  by  constructing  an  equivalent  model  based  on  the  beam  inputs  described 
above.  The  procedure,  illustrated  in  Fig.  4-2  for  Love  waves  into  a  downdip, 
is  to  replace  the  upstream  waveguide  (Fig.  4-2a)  by  the  equivalent  beam  input 
at  the  dipping  layer  junction  (Fig.  4-2b) .  By  extending  the  dipping  interface 
back  through  the  beam,  the  problem  is  analogous  to  a  single  beam  of  SR-rays 
projected  into  a  wedge  or  prism.  A  similar  picture  holds  for  Love  waves  into 
an  updip. 

In  this  context,  the  two-point  ray  tracing  problem  in  the  dipping  layer 
can  be  stated  as:  Given  a  receiver  point,  R  ,  in  the  prism,  find  the  inter¬ 
face  entry  points  of  those  rays  in  the  equivalent  beam  which  are  incident  on 
R  .  The  advantage  of  the  equivalent  model  is  that  the  prism  is  a  semi-infinite 
sector  and  lacks  a  characteristic  length.  The  only  lengths  involved  are  the 
beam  width  and  entry  point,  which  are  boundary  conditions  and  do  not  affect 
the  intrinsic  geometry.  The  virtue  of  this  situation  is  that  lacking  a  length 
scale  the  ray  system  must  exhibit  self-similarity,  simplifying  the  ray  picture 
considerably.  The  two-point  solution  follows  by  analyzing  the  rays. 
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4.1.1  The  self -similar  system  of  rays — Given  a  prism  with  included  angle, 

Y  ,  suppose  that  an  SH-ray  is  shot  in  the  downdip  direction  from  a  source  near 

the  vertex  at  radius  r^  and  takeoff  angle,  <5^  ,  on  the  lower  face.  Fig.  4-3. 

The  system  of  multiply-reflected  rays  is  immediately  drawn  on  the  basis  of 

straight  ray  paths  and  Snell's  law  for  internal  reflection  (reflection  angle 

equal  to  incidence  angle).  Reflections  occur  at  r  with  reflection  angle, 

n 

6  ,  n  ■  2,3,...  measured  from  the  normal  directed  into  the  prism.  Because  the 

n 

geometry  lacks  a  length  scale,  the  ray  system  must  exhibit  this  symmetry.  Dim¬ 
ensionally,  no  characteristic  length  requires  that  the  variables  and  parameters 
of  the  ray  system  reduce  to  r/r^  ,  9  ,  S  ,  iere  r  and  0  are  polar  coordi¬ 
nates.  Consequently,  a  solution  for  r/r^  the  ray  must  take  the  form 

r/t1  •  fCe^.y)  , 

i.e.  homogeneous  of  degree  one.  The  graphical  result  is  shown  in  Fig.  4-4  by 
extending  all  of  the  internal  rays  back  towards  the  vertex,  where  they  are 
clearly  tangent  to  a  single  circle  centered  at  the  vertex.  Radius  of  the  tan¬ 
gent  circle,  r^  depends  on  the  source  location  and  takeoff  angle  as 

r0  *  rlcos6l  *  <*-D 

A  little  geometry  shows  the  extended  rays  are  tangent  to  the  circle  at  the 
angles  shown,  hence 

Tq  m  r^co$6^  i  a  ■  1,2,  ...  ,  (4—2) 

a  very  powerful  result  for  the  two-point  problem.  The  formula  for  takeoff  angle 
in  (4-2)  is  found  from  the  geometrical  requirement  of  Snell's  law,  that  incidence 
angle  increases  by  the  prism  angle  from  one  reflection  to  the  next,  i.e.  the  re¬ 
cursion  6  -  6  +  Y  *  whence 

n+1  n 

5n  "  51  +  (n-1>Y  •  '  (4-3) 
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Clearly,  6q  muse  be  less  than  it/ 2  ,  which  places  an  upper  bound  on  the  num¬ 
ber  of  reflections.  Solving  6^  <  tt/2  for  n  yields 

it/2-6 

n  <  1  + 

from  which  the  upper  bound  is 


where  INT(  )  gives  the  integer  part  of  its  agreement.  The  Nth  reflection 
is  unique  in  that  the  reflected  ray  never  reaches  the  opposite  face,  but  goes 
to  infinity  into  the  prism  interior. 


4.1.2  Rays  to  an  interior  point — The  above  relations  characterize  the 

semi-infinite  ray  path  from  the  entry  point  to  reflection  points  to  the  point 

at  infinity  along  the  Nth  ray.  Next,  consider  the  finite  ray  path  from  the 

entry  point  to  some  interior  point  on  one  of  the  ray  segments  illustrated  in 

Fig.  4-5.  Assume  that  the  point  is  on  the  nth  crossing  ray  (n  even  or  odd) 

at  Rq  and  G  ,  where  the  angle  is  measured  from  the  upper  face  as  shown.  In 

the  figure,  applying  the  law  of  sines  to  the  triangles  connecting  the  vertex 

0  ,  r  on  a  face  and  R  in  the  interior,  radius  to  the  interior  point  is 
u  n 


(4-8) 


Length  of  the  ray  path  from  the  entry  point  to  the  reflection  at  r^  is 

L  ■  r  sin6  -  r,sin6.  , 
n  n  n  1  1 

which  follows  by  unfolding  the  ray  outside  the  prism  about  each  reflecting  face. 

Referring  back  to  Fig.  4-4,  these  unfolded  rays  are  coincident  with  the  extended 

ray  tangent  to  the  circle.  Images  of  the  corresponding  reflection  points  are 

A 

marked  r^  in  the  figure.  Therefore,  the  total  path  length  from  the  entry 

point  to  a  point  on  the  nth  crossing  ray  at  R  ,  Q  is 

i  n 

sin  8  "| 

S5^r+TT  '  risln6i  ’  (4“9) 

n  n  J 


L  +  l 


sin5  + 


where  r^  ,  and  0^  are  given  above.  Note  that  r^sin6^  is  just  the 
distance  of  the  entry  point  to  the  circle  tangency  point.  This  is  a  conveni¬ 
ent  point  of  constant  phase  for  the  equivalent  beam  input  described  earlier. 


4.1.3  The  two-point  solution— The  two-point  ray  tracing  problem  can  now 
be  solved.  The  solution  is  illustrated  in  Fig.  4-6.  Assume  that  the  receiver 
is  located  at  R  ,  0  and  that,  for  the  time  being,  single  ray  sources  are 
distributed  over  the  entire  lower  face.  Referring  to  the  figure,  there  are 
exactly  N  entry  points,  r^n  on  the  lower  face  which  illuminate  the  receiver. 
These  correspond  to  the  direct  ray,  n  ■  1  and  all  higher  multiples,  n  ■  2  , 

. . .  ,  N  ,  and  follow  by  rewriting  (4-5)  using  (4-2)  as 


R  cosiS. 

r7  “  f(Q;n)  S  cos(0  +6  ) 


'1  ~  '  n  n' 

Setting  R  ■  R  and  solving  for  the  corresponding  r. 

n  i 


(4-10) 

(renamed  r^R)  yields 


rln  “  f  (0;n)  ;  n  “  1,2 . N  *  (4-11) 

If  the  entry  point  is  within  the  equivalent  beam  input  window,  then  it 
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is  outside  the  win- 


correspoads  to  s  Love  ray  from  the  waveguide,  but  if  r^Q 
dow  it  has  no  significance,  at  least  in  the  present  context.  Therefore,  given 
an  input  window  corresponding  to  some  Love  wave,  all  rays  incident  on  the  re¬ 
ceiver  can  be  determined  from  <4-11) . 

4.1.4  Amplitude  and  phase  at  the  receiver — Receiver  response  is  the  sum 
of  the  incident  harmonic  waves  traveling  along  the  rays,  each  with  its  own  phase 
and  amplitude.  Consider  M  rays  illuminating  the  receiver,  each  with  amplitude 
A^  and  phase  P  ,  then  response  is 

M 

v(R,0,t)  -  l  A.  sin(uit-P.) 

i-1  1  1 


(4-12) 


(4-13) 


(4-14) 


The  individual  ray  amplitudes  and  phases,  and  are  calculated  on  the 

basis  of  total  ray  length  and  reflection  history.  The  amplitude  changes  at 
each  subcritlcal  reflection  (B^sin  80  that  A.  is  just  the  product 

of  the  entry  amplitude,  A°  ,  and  all  subcritlcal  reflection  coefficients,  R^  , 
j  ■  1  ,  ,  on  the  interface. 


A 


i 


(4-15) 
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The  phase  on  each  ray  is  determined  from  the  ray  length  and  phase  changes  at 
supercritical  interface  reflections. 


Ki  v 

l  *£  + 

1  k-1  X 


2ir  T 


(4-16) 


o  k, 

where  P ,  is  the  phase  at  the  entry  point,  ,  k»l  ,  K.,  are  the  phase  changes 

A 

at  supercritical  reflections,  is  the  total  ray  length  to  the  receiver,  and 

/s 

X  -  is  the  wavelength  in  the  prism.  The  are  measured  from  some 

convenient  point  of  constant  phase  in  the  equivalent  beam.  Reflection  coeffici¬ 
ents  may  be  found  in  Aki  and  Richards  [1980],  for  example. 


4.1.5  Updip  versus  downdip  ray  tracing' — Although  the  preceding  geometric 
ray  analysis  was  done  for  the  downdip  case,  it  is  equally  valid  for  updip.  This 
can  be  seen  by  continuing  the  entry  angle,  5^  ,  in  Fig.  4-3  from  positive  values 
(downdip)  through  zero  to  negative  values  (updip).  To  consider  the  updip  case, 
it  is  only  necessary  to  let  6^  be  negative  in  the  formulas,  by  analytic  con¬ 
tinuation.  This  yields  a  significantly  different  ray  picture  because  as  the  ray 
shoots  updip,  from  (4-3)  its  incidence  angle  increases  by  the  dip  angle  through 
negative  values  at  each  reflection  until  it  passes  through  zero.  At  this  point 
incidence  angle  is  positive  and  the  ray  shoots  downdip,  i.e.  the  ray  has  turned. 
Such  a  case  is  drawn  in  Fig.  4-7. 

In  application  to  Love  waves,  this  Implies  phenomenological  differences 
between  updipping  and  downdipping  transitions.  Rays  into  a  downdip  tend  to 
Increase  both  their  Incidence  angle  and  spacing  between  reflection  points. 

Fig.  4-3.  In  contrast,  rays  into  an  updip  decrease  their  incidence  angle  and 
reflection  spacing  up  to  the  turning  ray,  and  increase  them  thereafter  in  the 
downdip  direction,  Fig.  4-7.  As  a  consequence.  Love  rays  into  a  downdip  remain 
trapped  in  the  layer  by  total  reflection  because  the  angles  become  more  grazing; 


but  into  an  updip  the  steepening  rays  will  generally  reflect  partially  at  the 
Interface  as  the  turning  point  Is  approached,  thereby  leaking  much  of  their 
energy  Into  the  lower  half space. 


In  terms  of  the  equivalent  beam  input,  beam  segments  crossing  the  layer 
In  a  downdip  will  spread  apart  while  In  an  updip  they  will  overlap,  caused  by 
the  respective  increase  and  decrease  in  reflection  spacing  on  the  faces.  The 
spread  due  to  downdip  causes  shadows  to  appear,  i.e.  regions  not  illuminated  by 
first-order  rays  from  the  beam.  The  beam  is  trapped  and  after  the  Nth  reflec¬ 
tion  goes  to  infinity  into  the  interior.  These  effects  are  shown  in  Fig.  4-8 
for  a  7°  downdip  with  five  rays  across  the  incident  beam,  as  in  Fig.  4-2.  In 
contrast,  for  updip,  constructive  Interference  in  the  overlapping  regions  gives 
rise  to -amplification,  and  when  the  beams  reflect  partially  off  the  interface 
they  transmit  into  the  halfspace.  Bay  beams  into  a  3s  updip  with  B^/B^  ■  2 
are  drawn  in  Fig.  4-9.  All  leaked  beams  are  Included  although  only  the  first 
few  have  significant  amplitude. 

Finally,  note  that  beams  can  propagate  no  further  into  the  transition  than 
the  turning  point  of  their  rightmost  ray  for  updip;  while  for  downdip  the  beam 
cannot  reach  the  free  surface  past  the  N  or  N-l  reflection.  Phase  velocities 
on  the  free  surface  decrease  down  the  transition,  approaching  B^  for  downdip, 
but  increase  for  updip,  becoming  high  near  the  turning  point. 

4.2  Results  of  Love  Wave  Propagation  in  Sedimentary  Prisms 

Response  at  a  receiver  in  the  prism  can  now  be  determined  for  Love  waves 
incident  from  the  uniform  waveguide.  First  the  dispersion  relation  for  a  par¬ 
ticular  model  and  mode  is  solved  for  a  discrete  set  of  phase  velocity-period 
(V-T)  pairs,  which  are  converted  to  ray  angle-period  (9L~T)  pairs.  Then  from 
the  equivalent  ray  beam  input  (Fig.  4-2),  response  to  first-order  rays  is  found 
for  each  0^-T  pair  using  the  recursions  and  sums  derived  above  for  two-point 


ray  tracing  in  the  prism. 
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Two  Love  wave  models  will  be  investigated  representing  sedimentary  wave¬ 
guides  2  km  and  5  km  thick  on  the  eastern  and  western  flanks  of  the  Rocky  Moun¬ 
tains,  respectively.  Approximate  attenuation  is  Included  to  indicate  the  effect 
of  material  damping  on  response.  Three  modes  will  be  examined  for  each  model. 
The  transitions  include  updipplng  and  downdipping  layers  at  3°,  5°  and  10°. 
Response  in  the  transition  is  evaluated  at  an  array  of  discrete  points,  either 
over  the  free  surface,  X  or  through  the  depth,  Z  .  For  each  Love  wave  mode  in 
a  model,  displacement  amplitude  will  be  calculated  both  as  a  function  of  Love 
wave  period  and  distance  or  depth.  Results  are  presented  in  the  form  of  three- 

A 

dimensional  surface  plots  typically  involving  10  discrete  amplitudes  per  plot 
(100  receivers  x  100  periods)  to  achieve  adequate  resolution. 

4.2.1  Ideal  free  surface  response — A  typical  example  of  surface  response 
is  shown  in  Fig.  4-10  for  the  2  km  waveguide  into  5”  updipplng  and  downdipping 
transitions.  Amplitude  is  a  function  of  Love  wave  period  and  distance  from  the 
end  of  the  uniform  waveguide  into  the  transition,  and  is  normalized  by  amplitude 
in  the  waveguide.  In  addition  to  the  response  surface,  peak  amplitude  over  the 
transition  as  a  function  of  period  is  plotted  in  the  vertical  amplitude-period 
plane.  The  plots  graphically  illustrate  the  fundamental  difference  between  up- 
dlp  and  downdip  propagation  in  a  transition.  Namely,  updip  is  characterized  by 
patterns  of  amplification  and  downdip  by  patterns  of  deampliflcation.  The  ex¬ 
planation  follows  from  the  geometric  comparison  of  updip  and  downdip  in  Sec. 
4.1.5. 

Consider  updip  response  for  a  particular  period  in  Fig.  4-10.  Moving  along 
the  transition,  near  the  origin,  response  is  identical  to  that  in  the  waveguide 
until  the  point  of  first  beam  overlap  is  reached,  where  the  amplitude  doubles. 
Past  the  overlap  the  amplitude  is  unity  until  the  next  sequence  of  beam  overlaps 
amplifies  it  by  constructive  interference.  Continuing  past  the  point  of  maximum 
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amplification,  amplitude  decreases  roughly  stepwise  and  then  very  quickly 
to  zero  as  the  bi,  •  is  turned  in  the  downdip  direction.  Interface  reflections 
near  the  turning  point  are  partial,  transmitting  most  of  the  incident  energy 
into  the  half space.  Consequently,  backscattered  rays  have  negligible  amplitude 
and  effect. 

For  downdip  response,  again  consider  a  particular  period.  Near  the  origin, 
response  is  identical  to  that  in  the  waveguide  up  to  the  first  shadow  where  the 
amplitude  drops  to  zero.  This  is  followed  by  a  zone  of  unit  response,  shadow, 
etc.  These  alternations  are  caused  by  gaps  between  the  ray  beam  on  the  free 
surface  as  it  reverberates  in  the  downdip  direction. 

The  patterns  of  abrupt  variation  in  amplitude  for  both  updip  and  downdip 
are  caused  by  either  a  change  in  the  number  of  rays  (beams)  incident  on  the 
receiver  point;  or  a  change  in  one  or  more  interface  reflections  on  a  ray,  from 
partial  to  total  or  vice  versa.  In  both  dip  cases,  the  sequence  of  jumps  is 
seen  to  smoothly  migrate  down  the  transition  as  the  period  decreases  due  to  the 
corresponding  increase  in  Love  ray  angle. 

To  track  the  amplitude  discontinuities  in  the  X-T  plane,  it  is  convenient 
to  examine  the  downdip  case.  By  simple  geometry  the  first  shadow  demarcation 
appears  at  distance,  X^  ■  h  tan  @L(T)  from  the  waveguide.  Note  that  this  is 
half  the  equivalent  input  beam  width.  The  next  demarcation  can  likewise  be  found 
at  X2  *  h  tan  (0^(T)  +  2y) ,  etc.  As  9^  becomes  larger  (grazing  rays)  for 
shorter  periods,  X^  and  X^  increase  quickly  as  the  argument  of  the  tan  func¬ 
tion  approaches  ir/2  .  This  explains  the  smooth  sweep  in  Fig.  4-10  nearly  para¬ 
lleling  the  T  axis  for  longer  periods  and  then  paralleling  the  X  axis  for 
shorter  periods. 

4.2.2  A  comparison  of  response  for  both  models — Plots  of  surface  response 
over  both  sedimentary  models  are  shown  in  Figs.  l-8a,b  for  four  cases  of  updip 
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and  downdip.  In  most  part  they  exhibit  the  same  behavior  as  described  above 
for  Fig.  4-10.  The  resemblance  arises  from  the  essential  similarity  of  the 
fundamental  mode  dispersion  curves.  Peak  updip  amplifications  are  comparable 
for  the  tvo  models  because  of  similar  Impedance  contrasts  between  half space  and 
layer  measured  by  their  ratio,  ■  1.88  and  2.2  for  the  5  km  and  2  km 

layers.  There  is  virtually  no  backscattering  of  energy  from  the  updips  because 
of  these  low  impedance  contrasts.  However,  greater  contrasts  will  produce  sig¬ 
nificant  ref lectiohs  back  into  the  waveguide.  Comparing  the  3°,  5°  and  10°  cases 
shows  that  the  discontinuous  surface  response  becomes  sparser  as  dip  angle  in¬ 
creases,  as  well  as  decreasing  peak  amplitudes.  The  reason  is  fewer  beam  rever¬ 
berations  at  the  greater  dip  angles. 

4.2.3  Higher  modes  and  response  through  depth — The  fundamental  mode  and 
its  first  two  overtones  were  calculated  for  both  models.  Their  features  are 
illustrated  in  Fig.  4-11  for  the  5  km  layer  with  5°  updip  and  downdip.  The 
fundamental  mode  exists  over  all  periods  while  the  overtones  exhibit  cutoffs 
which  are  evident  in  the  picture.  Correlation  of  overtone  response  with  that 
of  the  fundamental  can  be  made  using  the  dispersion  curves.  These  show  that 
shape  of  the  response  surface  is  a  function  principally  of  the  Love  wave  ray 
angle  (i.e.  phase  velocity).  In  fact,  a  stretching  transformation  on  overtone 
period,  to  allign  the  overtone  dispersion  curves  with  the  fundamental  mode 
curve,  provides  the  mapping  between  response  curves.  This  is  evident  in  a 
comparison  of  mode  2  and  mode  3  response,  e.g.  stretch  the  T  scale  to  double 
its  present  length  on  the  third  mode  and  compare  to  the  second. 

Mode  shapes  for  the  fundamental  and  two  overtones  are  shown  in  Fig.  4-12 
for  the  5  km  layer.  These  are  merely  plots  of  amplitude  versus  depth  and  period 
at  the  waveguide  junction  into  an  arbitrary  downdip.  They  were  compared  to  the 
analytic  expression,  Icosoiq^zj  to  verify  the  analysis.  Amplitude  at  the 
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interface  depth  in  the  figure  demonstrates  the  relative  strength  of  the  evanes¬ 
cent  wave  in  the  halfspace.  Clearly,  the  evanescent  waves  are  stronger  as  the 
interface  amplitude  increases  at  the  longer  periods. 

Amplitude  versus  depth  in  the  transition  is  plotted  in  Fi».  4-13  for  2  km, 
5s  updip  and  downdlp.  The  corresponding  surface  response  was  described  in 
Sec.  4.2.1.  Downdip  cross  sections  are  at  10  and  20  km  from  the  waveguide,  and 
updip  sections  are  at  5  and  10  km.  These  plots  show  that  for  downdlp  the  rever¬ 
berating  but  non-overlapping  beam  creates  a  diamond  pattern  of  steps  in  the 
Z-T  plane.  The  amplitude  is  zero  in  the  shadows,  .5  where  one  ray  is  incident 
and  1.  where  two  rays  intersect.  Spacing  along  the  period  axis  decreases  with 
distance  from  the  waveguide.  For  updip,  the  diamond  patterns  are  still  evident 
but  amplification  due  to  overlapping  complicates  the  picture.  It  appears  that 
amplification  on  the  surface  implies  the  same  through  depth,  although  the  ampli¬ 
tude  is  not  a  smooth  function  of  depth  as  for  the  incident  mode.  Fig.  4-12. 
Higher  modes  exhibit  similar  patterns. 

4.2.4  Effects  of  material  damping  on  surface  response — The  analysis  des¬ 
cribed  in  Section  4.1  was  for  waves  in  ideal,  non-dissipative  geologic  media. 

In  reality  such  materials  exhibit  material  damping.  To  indicate  what  effects 
this  damping  might  have  on  transition  response,  a  frequency  independent  Q  was 
introduced  in  the  ray  analysis  by  merely  attenuating  amplitude  on  the  basis  of 
2it/Q  <■  -2AA/A  ■  -2d  ,  where  AA  is  the  change  in  amplitude  for  one  cycle  on 
the  ray  and  6  is  the  logarithmic  decrement.  This  is,  of  course,  a  crude 
approximation  but  for  low  damping  values  it  is  not  unreasonable,  see  Aki  and 
Richards  [1980]. 

Both  models  were  examined  for  a  number  of  damping  values.  The  results 
already  shown  were  for  Q  -  1000  ,  essentially  undamped.  Effects  are  shown 
in  Fig.  4-14  for  the  2  km  shale /sandstone  waveguide,  with  Q  ■  1000,100,40  . 
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Clearly,  high  frequencies  are  most  highly  attenuated  but  response  Is  still  signi¬ 
ficant  at  moderate  distances  (10-20  waveguide  depths). 


4.3  Extension  of  Ray  Tracing  to  P-SV-Waves  in  a  Sedimentary  Prism 

The  analysis  of  Rayleigh  (P-SV)  waves  in  a  sedimentary  prism  is  based  on  an 
extension  of  the  theory  presented  in  Section  4.1.  The  Rayleigh  case  involves 
two  wave  types,  P-  and  SV-body  waves,  whereas  the  Love  case  involves  only  SH- 
waves.  This  section  discusses  two  aspects  of  the  Rayleigh  wave  case  requiring 

f 

special  attention.  First,  an  algorithm  is  needed  to  deal  with  the  large  number 
of  internal  reflections  generated  in  a  sedimentary  prism.  Second,  a  method  is 
needed  to  track  the  energy  which  leaks  into  the  underlying  halfspace  due  to 
subcritical  interface  reflections;  this  has  practical  application  to  tracking 
surface  wave  energy  leaked  from  basin  edges. 

4.3.1  Binary  tree  representation  of  rays — Analysis  of  the  Rayleigh  case 
requires  more  elaborate  bookkeeping  because  at  each  reflection  both  P-  and  SV- 
waves  are  usually  generated;  in  fact,  the  number  of  rays  multiply  exponentially 
and  quickly  overwhelm  the  storage  capacity  even  of  large  computers.  This  pro¬ 
blem  is  solved  by  screening  out  low  amplitude  rays. 

A  general  scheme  for  calculating  all  P-SV-rays  is  required.  The  scheme 
described  here  involves  recursive  application  of  the  canonical  solution  for 
ray  tracing  in  a  prism.  Figure  4-15  illustrates  a  typical  P-SV-ray  system  in 
a  layer.  The  system  consists  of  a  main  trunk  (the  incident  ray)  from  which  a 
cascade  of  branches  (the  reflected  rays)  are  produced.  Each  branch  in  turn 
may  be  regarded  as  a  trunk.  This  situation  is  conveniently  represented  by  a 
binary  tree.  The  term  "binary"  refers  to  the  fact  that  a  single  incident  ray 
usually  produces  two  rays  upon  reflection.  The  binary  tree  data  structure  is 
especially  convenient  for  keeping  account  of  ray  data  (amplitude,  phase)  in  a 
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format  that  makes  screening  easy  and  computationally  economical.  In  this  case, 
we  screen  out  rays  whose  amplitudes  are  below  a  threshold  value. 

To  illustrate  how  the  binary  ray  tree  is  used  to  screen  rays.  Fig.  4-lGa 
depicts  a  fan  of  rays  emanating  from  a  point  on  a  dipping  interface.  Some  of 
the  rays  (positive  angle)  are  in  the  downdip  direction  while  others  (negative 
angle)  are  in  the  updip  direction.  The  binary  tree  for  each  ray  is  computed 
as  follows.  All  reflection  points  and  reflection  coefficients  are  calculated 
for  the  trunk  ray.  When  a  branch  occurs  (e.g.  when  a  P-wave  is  reflected  as 
S-  and  P-wa/es) ,  its  amplitude  is  calculated  and  compared  to  the  threshold  value. 
If  the  branch  amplitude  is  significant  (i.e.  greater  than  the  threshold)  the 
branch  data  is  retained  in  a  memory  stack;  if  the  branch  amplitude  is  below  the 
threshold,  the  branch  is  pruned  and  no  further  calculations  using  it  are  made. 
This  sequence  is  followed  down  the  main  trunk  until  all  branches  are  accumulated 
or  discarded.  The  first  retained  branch  off  the  trunk  then  is  treated  as  if  it 
were  a  trunk  and  procedure  is  repeated.  The  second  and  subsequent  branches  are 
treated  similarly. 

A  graphical  illustration  of  binary  tree  data  in  which  P-  and  S-rays  are 
shot  at  various  updip  and  downdip  angles  is  given  in  Figs.  4-17  and  4-18.  No 
pruning  is  applied  to  these  trees.  As  defined  in  Fig.  4-16b,  S-rays  are  drawn 
to  the  left  and  P-rays  to  the  right.  P-branches  generally  undergo  more  reflec¬ 
tions  than  S-branches.  In  Fig.  4-17,  for  example,  each  reflection  is  denoted 
by  a  unit  of  vertical  growth  of  the  branch.  Clearly,  branches  tend  to  grow 
more  to  the  right  than  to  the  left.  Ray  trees  for  S-rays  are  shown  in  Fig.  4-18. 
In  cases  of  shallow  updip,  the  initial  reflections  of  the  S-trunk  are  so  shallow 
that  no  P-branches  are  initially  generated.  For  steeper  incidence  cases,  both 
P-  and  S-branches  are  generated  at  the  first  reflection. 
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4.3.2  Leakage  of  surface  wave  energy  at  basin  edges — Subcritical  reflec¬ 


tions  at  the  interface  between  the  sedimentary  prism  and  underlying  bedrock 
result  in  leakage  of  surface  wave  energy  out  of  the  prism.  This  is  illustrated 
in  Fig.  1-9  for  Love  waves  of  three  different  periods  propagating  into  a  3°  up- 
dipping  transition.  In  terms  of  our  equivalent  ray  beam  analysis,  the  energy 
is  leaked  as  ray  beams  at  angles  depending  on  the  dip  angle,  impedance  contrast 
and  period.  This  treatment  neglects  diffraction  effects  which  would  cause  beam 
spreading  and  would  modify  amplitudes.  To  include  diffraction  effects  requires 
numerical  analysis,  which  is  the  subject  of  current  research. 

A  similar  mechanism  occurs  in  the  case  of  Rayleigh  wave  overtones  propagating 
into  an  updip;  however,  in  this  case,  both  P-  and  SV-beams  are  leaked.  An  example 
for  a  single  ray  incident  on  a  10°  updip  is  shown  in  Fig.  4-19.  The  leaked  rays 
are  shown  in  Fig.  4-t9a  for  an  incident  P-ray  at  three  different  angles.  Two 
numbers  describe  each  ray — the  upper  is  the  amplitude  and  the  lower  is  the  angle 
of  refraction  measured  from  the  vertical.  An  important  finding  of  the  analysis 
Is  that  ray  beams  of  significant  amplitude  leave  the  prism  at  a  variety  of  angles. 
As  a  result,  these  beams  will  re-emerge  at  the  surface  over  a  wide  range  of  dis¬ 
tances  from  the  prism.  This  is  illustrated  in  Fig.  4-2D  where  a  fan  of  rays, 
each  of  which  may  correspond  to  a  beam  leaked  at  an  angle  between  34s  and  55° 
into  the  Sinclair,  IVY  crust,  re-emerges  at  widely  varying  distances  from  the 
point  of  incidence.  Rays  leaked  at  angles  less  than  34s  or  more  than  55*  would 
re-emerge  at  distances  greater  than  260  km  or  less  than  10  km,  respectively.  The 
leaked  rays  for  at-,  incident  SV-ray  are  shown  in  Fig.  4-19b,  where  the  findings 
with  respect  to  the  variety  of  exit  angles  are  similar.  Comparing  Fig.  4-19a 
and  b  shows  that,  in  general,  S-rays  will  propagate  further  up  the  transition 
than  P-rays.  Further  study  of  this  problem,  including  diffraction  effects,  is 
warranted . 
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Figure  4-2.  The  equivalent  Love  wave  transition  model 
based  on  a  ray  beam  projected  Into  a  prism 
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The  geometry  of  rays  incident  on  an 
interior  point  at  ,  0  . 
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Figure  4-7. 


The  complete  ray  system  for  a  ray  shot 
in  the  updip  direction. 
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Figure  4-8.  The  equivalent  Love  ray  beam 
in  a  7*  dovndip. 
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Figure  4-10.  An  example  of  the  2  km  waveguide  fundamental  Love  mode 
surface  response  as  a  function  of  wave  period  and 
distance  over  the  transition. 
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Figure  4-14.  Effects  of  material  damping  on  2  km,  3°  updip  and  downdip 
surface  response. 
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Figure  4-15.  A  typical  P-SV-ray  cascade  in  a  layer 
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.  A  ray  fan  in  the  prism 


b.  The  P-  and  S-ray  convention  in  the  binary  tree 


Figure  4-16.  The  ray  fan  in  a  prism,  and  the  convention  for 
drawing  the  binary  tree  for  each  ray  of  the  fan. 


f 


It1'1 


Incident 
P-ray  — ' 


(M8  Y0.14\ 
-14.2  -24.4 


Figure  4-1 9a. 


Leaked  raya  from  a  10°  prism  for  an  incident  P-ray. 
Upper  number  is  amplitude  and  lower  is  angle  from 
vertical  in  degrees  (negative  sign  for  updip).  Rays 
In  prism  are  pruned  to  an  amplitude  of  10X  of  inci¬ 
dent  ray. 


SECTION  5 


SUMMARY  AND  CONCLUSIONS 

This  section  briefly  summarizes  the  work  performed  under  the  present 
contract  and  presents  the  main  conclusions  prompted  by  three  questions  stated 
in  Section  1.5.  The  methods  used  in  this  study  consider  Idealization  of  the 
geology  and  a  portion  of  the  complete  wave  field.  To  consider  more  complete 
solutions  requires  full  wave  field  numerical  analysis. 

5.1  Summary 

This  report  examines  the  propagation  of  surface  waves  from  the  Pocatello 
Valley,  Idaho  earthquake  sequence  across  Wyoming  and  the  central  Rocky  Noun* 
tains  to  Wing  V.  A  model  of  the  travel  path  between  Pocatello  Valley  and  Wing  V 
is  prepared  on  the  basis  of  geological  reconnaissance  including  geology  and  tec¬ 
tonics  of  the  Rocky  Mountain  region.  Cross  sections  along  the  travel  path  show 
basin  sediments  overlying  basement  rock  thinning  toward  the  east  with  interface 
dips  of  2a-10a.  A  range  of  crustal  velocity  models  were  found,  including  Flam¬ 
ing  Gorge,  Sinclair  and  Wolcott,  and  were  augmented  with  data  from  oil  well 
sonic  logs  for  near-surface  sediments.  The  arrangement  of  principal  uplifts 
and  basins  in  the  region  is  three-dimensional  and  cannot  be  idealized  by  simple 
one-  or  two-dimensional  travel  path  models. 

Short-  and  long-period  seismograms  recorded  at  Golden,  Colorado  during  the 
1975  Pocatello  Valley  earthquake  were  examined  to  infer  bedrock  motion  near 
Wing  V.  The  main  event  was  not  well  recorded  but  a  magnitude  4.7  aftershock 
provided  data  on  body  end  fundamental  mode  Rayleigh  wave  phases.  From  the  long- 
period  surface  wave  data,  a  group  velocity  dispersion  curve  and  particle  motion 
orbits  (elllpticlty)  were  obtained.  Theoretical  dispersion  curves  in  crustal 
models  along  the  travel  path  were  compared  with  the  one  derived  from  the  Golden 


records.  The  comparison  shows  that  the  measured  dispersion  curve  for  the  range 
of  periods  5-12  sec  lies  within  the  bounds  obtained  theoretically  using  the 
various  models  along  the  travel  path.  Favorable  agreement  between  measured 
and  calculated  elliptlcity  of  surface  displacements  for  fundamental  mode  Ray¬ 
leigh  wave  was  also  obtained.  This  gives  us  confidence  in  extrapolating  the 
theoretical  model  into  the  period  range  of  1-5  sec,  which  is  significant  for 
Wing  V  seismic  alarms.  In  all  theoretical  cases  examined,  group  velocity 
minima  were  found  in  this  shorter  period  range  suggesting  that  Airy  phases 
capable  of  producing  long  duration  ringing  of  the  crust  may  be  significant. 

There  are  two  classes  of  surface  waves  which  are  incident  at  Wing  V. 

One  is  the  fundamental  Rayleigh  or  free  surface  wave  consisting  of  interfering 
evanescent  waves.  The  other  is  the  channeled  surface  wave  including  Rayleigh 
overtones  and  Love  wave  modes,  all  of  which  are  composed  fully  or  in  part  of 
trapped  plane  waves.  This  channeled  type  may  be  readily  analyzed  using  con¬ 
ventional  ray  theory  in  contrast  to  free  surface  waves.  In  preparation  for 
the  study  of  channeled  wave  Interaction  with  basin  and  mountain  structure,  we 
examine  channeled  surf  ace  waves  using  seismic  ray  theory  in  smooth  crustal 
models.  The  results  show  ray  path  penentration  depths  and  travel  times  for 
P-rays.  The  finding  is  that  the  first  Rayleigh  overtone  (second  mode)  and 
the  first  two  Love  wave  modes  are  propagated  near  the  free  surface,  hence 
interact  with  lateral  variations  in  the  sedimentary  cover.  However,  higher 
modes  travel  deeper  in  the  crust  and  interact  much  less  with  surface  sedi¬ 
ments.  This  ray  viewpoint  also  provides  a  basis  for  representing  higher 
% 

mode  surface  waves  by  ray  beams  in  the  next  part  of  the  study. 

To  examine  the  behavior  of  channeled  surface  waves  along  the  travel 
path,  the  gently  dipping  sedimentary  layer  on  the  western  flank,  of  the  Rockies 
is  Idealized  as  a  homogeneous  prism  and  analyzed  by  a  two-point  ray  tracing 


scheme.  Input  to  the  prism  consists  o£  various  Love  modes  decomposed  into 
SH-ray  beams  mentioned  above  and  the  solution  is  constructed  by  tracking 
all  rays  Incident  on  an  array  of  receivers  at  varying  distances  along  the 
prism.  Examples  considered  are  sedimentary  prisms  of  5  km  and  2  km  depth 
at  the  input  point  for  various  updips  and  dovndips;  these  are  inspired  by 
geology  on  the  eastern  and  western  flanks  of  the  Rocky  Mountains ,  respectively. 
The  model  was  extended  to  consider  the  effects  of  typical  damping  values  for 
basin  sediments.  Overall,  these  results  show  that  little  surface  wave  energy 
is  reflected  at  the  flank  but  instead  is  transmitted  into  the  crust  at  a  vari¬ 
ety  of  angles.  To  include  Rayleigh  overtones  requires  ray  tracing  for  both  P- 
and  SV-rays  which  is  much  more  complicated  than  the  SH-wave  case  considered 
above.  As  a  step  toward  accomplishing  this,  a  binary  tree  representation  of 
P-SV-ray  tracing  was  developed  but  the  complete  ray  tracing  scheme  has  not  yet 
been  implemented.  Some  examples  showing  the  leakage  of  P-SV-rays  from  the  prism 
are  analyzed. 

5.2  Conclusions 

Here  we  describe  conclusions  reached  regarding  the  three  questions  raised 
in  Section  1.5.  The  first  question  concerns  wave  types  causing  seismic  alarms 
at  Wing  V.  For  any  surface  structure  subject  to  seismic  ground  shaking,  three 
important  factors  governing  response  are  amplitude,  frequency  and  duration  of 
shaking.  Because  the  Wing  V  missile  suspension  system  has  three  lightly  damped 
(less  than  l.Z  critical  damping)  resonances  in  the  period  range  1-5  sec,  input  of 
long  duration  ground  shaking  at  appropriate  periods  will  cause  seismic  alarms 
even  though  input  amplitudes  are  low.  Thus,  the  important  factors  are  proper 
frequency  and  sufficiently  long  duration;  for  example,  periods  of  about  3.5  sec 
and  durations  on  the  order  of  2  minutes. 
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Short  period  Golden  seismograms  from  the  1975  Pocatello  Valley  earthquake 


sequence  indicate  the  P  ,  P  and  P  are  of  too  short  duration  to  be  the 

8 

principal  cause  of  alarms.  However,  later  arrivals,  including  S  body  waves  . 
and  1.^  Love  and  Rayleigh  surface  waves,  possess  periods  and  sufficient 
durations  to  cause  alarms.  This  is  in  contrast  to  the  1979  St.  Elias,  Alaska 


earthquake  which  also  caused  seismic  alarms  due  to  first  arriving  P-waves  with 
periods  on  the  order  of  3  sec,  duration  on  the  order  of  five  minutes  and  peak 
amplitude  near  0.1  mm. 

The  travel  path  between  Pocatello  Valley  and  the  Denver  Basin  includes 
350  km  of  the  Green  River  Basin  with  5-7  km  of  sedimentary  cover  (half  the 
travel  path).  For  the  period  range  1-5  sec,  first  and  second  mode  Rayleigh 
waves  and  first  mode  Love  waves  are  confined  to  this  sedimentary  layer.  In 
contrast,  higher  mode  surface  waves  propagate  much  deeper  into  the  crust. 
Whereas  the  sedimentary  surface  waves  continually  interact  with  near-surface 
inhomogeneities,  the  higher  mode  crustal  surface  waves  are  affected  much  less 
because  most  of  their  component  ray  path  is  at  depth.  Our  ray  analysis  shows 
how  fundamental  mode  Love  waves  and  first  overtone  Rayleigh  waves  Interact  with 
the  sedimentary  prism  on  the  western  flank  of  the  Rocky  Mountains.  The  result 
of  this  interaction  is  leakage  of  energy  in  the  form  of  deeper  penetrating 
waves  which  may  or  may  not  re-emerge  as  body  or  higher  mode  surface  waves  near 
Wing  V.  The  higher  mode  crustal  surface  wave  phases  which  are  relatively  un¬ 
affected  by  the  travel  path  are  steeply  incident  at  Wing  V  and  in  terms  of 
their  ray  components  virtually  are  indistinguishable  from  body  waves  except 
for  time  of  arrival.  Therefore,  the  higher  mode  crustal  phases  in  the  vicinity 
of  Wing  V  can  be  viewed  as  steeply  Incident  body  waves  which,  as  shown  in  our 
previous  report,  Wojclk  et  al.  [1980],  are  responsible  for  anomalous  patterns 
of  seismic  alarms  due  to  resonance  of  the  underlying  sedimentary  column. 
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The  second  question  deals  with  the  response  anomalies  due  to  propagation 
of  surface  waves  through  sedimentary  prisms.  We  find  response  of  the  type 
illustrated  in  Fig.  4-10.  These  figures  show  amplifications  and  distinct 
cutoff  distances  for  propagation  into  an  updip,  and  nominal  surface  wave  re¬ 
sponse  interrupted  by  shadow  zones  for  propagation  into  a  downdip.  These  zones 
propagate  with  an  unsteady  speed  depending  on  the  dispersive  character  of  the 
Incident  surface  wave;  by  dispersive  character  we  mean  chiefly  the  rate  at 
which  surface  wave  period  changes  at  the  receiver  due  to  dispersion  of  the 
source  pulse.  These  results  are  relevant  to  siting  of  facilities  at  basin 
edges.  Due  to  the  distinct  distance  cutoff  for  updip  propagation,  facilities 
can  be  shielded  from  longer  period  components  of  surface  waves.  Although 
shorter  period  components  are  incident  throughout  the  prism,  they  tend  to  be 
highly  damped  in  sedimentary  and  basin  fills. 

The  third  question  is  concerned  with  leakage  of  surface  wave  energy.  In 
basins  where  depth  changes  are  continuous  with  interface  slopes  less  than  15° 
or  so,  we  find  that  propagation  over  a  downdip  tends  to  promote  shallowing  of 
the  component  rays  and  energy  trapping.  In  contrast,  propagation  over  an  up¬ 
dip  tends  to  promote  steepening  of  the  rays,  some  or  all  of  the  energy  leaked 
depending  on  the  geometry,  period  and  interface  Impedance  contrast.  The  leaked 
energy  can  be  expressed  as  spreading  ray  beams  projected  into  the  crust,  which 
will  re-emerge  at  some  distance  depending  on  the  refraction  angle  and  details 
of  the  velocity  variation  and  depth.  This  leaked  energy  may  have  the  same 
depth  of  penetration  (ray  parameter)  as  higher  mode  surface  waves  which  did 
not  Interact  with  the  prism  and  thereby  reinforce  them;  or  the  energy  may  re- 
emerge  near  the  receiver  as  a  spreading  beam  of  body  waves.  The  former  inter¬ 
pretation  is  appropriate  to  longer  distances  (multiple  skips)  between  source 
and  receiver;  the  latter  is  more  appropriate  at  shorter  distances  (single  skip). 
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